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PROCEDES ET DISPOSITIFS DE FILTRAGE NUMERIQUE 

La presente invention se rapporte a des procedes et a des dispositifs 
de filtrage numerique. 
5 Elle trouve une application avantageuse au filtrage d'images 

numeriques, en particulier pour la compression d'images conformement a 
"JPEG 2000", en cours de normalisation. 

On connaTt notamment deux techniques pour realiser (en anglais 
''implement) le filtrage de signaux numeriques en utilisant des fiitres en 
10 ondelettes (en anglais "wavelet filters") : 

- le filtrage par convolution et 

- le filtrage par augmentation ou "lifting". 

On rappelle ci-apres quelques notions elementaires sur les fiitres en 
ondelettes, puis sur le lifting. 
15 On sait que le filtrage en ondelettes d'un signal monodimensionnel Xj 

peut etre defini comme : 

- la convolution d'un filtre passe-bas H 0 (k) avec le signal Xj, le signal 
Y n resultant de la convolution etant sous-echantillonne (ou decime, en anglais 
"downsampled*) par un facteur de deux : 

20 Y 2n =£H 0 (k).X 2n _ k (1) 

k 

- la convolution d'un filtre passe-haut H^k) avec le signal Xj, le 
signal Y n resultant de la convolution etant decime par un facteur de deux : 

Y 2n+1 = X H i( k )- X 2n+1-k (2) 
k 

25 Le resultat du filtrage, appele operation de filtrage en ondelettes (Ho, 

Hi), produit un signal Yj dont chaque echantillon de rang pair Y 2n est un 
coefficient passe-bas et chaque echantillon de rang impair Y 2n +i est un 
coefficient passe-haut. 

Le filtre passe-bas et le filtre passe-haut doivent satisfaire certaines 

30 conditions, appelees conditions de reconstruction parfaite, de fa<?on que la 
reconstruction parfaite du signal X, a partir du signal Yj soit possible. 
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Si on utilise la notation de la transformee en z pour les filtres a 
reponse impulsionnelle finie, a savoir, 

Hj(z) = X H i( k )- z ~ k ' (3) 
k 

alors la condition de reconstruction parfaite que les filtres doivent satisfaire est : 

5 H 0 (z).H 1 (-z- 1 )-fH 0 (-z).H 1 (z- 1 ) = 2. (4) 

On sait que certains filtres en ondelettes presentent des proprietes 
de symetrie. 

Un filtre en ondelettes symetrique de longueur impaire H(k) (en 
anglais "odd-length symmetric wavelet fitted) ou filtre WSS (en anglais "Whole- 
10 Sample-Symmetric") est defini par la relation suivante, dite symetrie WSS, entre 
tous les coefficients du filtre : 

Vk, H(-k) = H(k) (5) 
Un filtre en ondelettes symetrique de longueur paire (en anglais 
"even-length symmetric wavelet filter) ou filtre HSS (en anglais "Half-Sample- 
15 Symmetric") est defini par la relation suivante, dite symetrie HSS, entre tous les 
coefficients du filtre : 

Vk, H(-k) = H(k + 1). (6) 
Un filtre en ondelettes antisymetrique de longueur paire (en anglais 
" even -l eng th antisymmetric wavelet filter 11 ) ou filtre HSA (en anglais "Half- 
20 Sample-Antisymmetric") est defini par la relation suivante, dite symetrie HSA, 
entre tous les coefficients du filtre : 

Vk, H(-k) = -H(k + 1). (7) 
La technique du lifting, mentionnee plus haut, consiste a realiser (en 
anglais "implement 1 ) une etape de filtrage en la decomposant en une sequence 
25 equivalente d'operations elementaires de filtrage, appelees etapes de lifting (en 
anglais "lifting steps"). Cette decomposition est appelee factorisation lifting (en 
anglais "lifting factorization"). Pour plus de details sur la factorisation lifting, on 
se reportera utilement a I'article de Ingrid DAUBECHIES et Wim SWELDENS 
intitule "Factoring wavelet transforms into lifting steps", in J. Fourier Anal. Appl., 
30 vol. 4, n° 3, pages 247 a 269, 1998. 
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On rappelle en particulier que, lorsqu'on realise une operation de 
filtrage sous forme d'etapes de lifting, pour transformer un signal Xj en un signal 
Yj, on commence generalement par copier tous les echantillons du signal Xj 
dans le signal Yj (i.e. Yj = Xj), puis on modifie alternativement les echantillons 
5 de rang pair et les echantillons de rang impair du signal Yj. 

II existe deux types d'etapes de lifting : 

- les etapes de lifting passe-haut et 

- les etapes de lifting passe-bas. 

Une etape de lifting passe-haut consiste a modifier des echantillons 
10 de rang impair (ou d'indice impair, en anglais "odd-indexed samples") (passe- 
haut) en leur ajoutant une fonction d'echantillons de rang pair (ou d'indice pair, 
en anglais "even-indexed samples") ponderes ou filtres : 

r \ 



Vn ' Y 2n + 1 = Y 2n+1 + R 



X a k Y 2n+2k 

V k J 



(8) 



La fonction R est, de fa?on generale, toute approximation de la 
15 variable x. La fonction R peut etre un operateur d'arrondi qui arrondit une valeur 
reelle x a un entier (tel que rentier le plus proche), ou peut etre simplement 
I'identite : R(x) = x. 

Une etape de lifting passe-bas consiste a modifier des echantillons 
de rang pair (passe-bas) en leur ajoutant une fonction d'echantillons de rang 
20 impair ponderes ou filtres : 



Vn,Y 2n =Y 2n +R 



X b k- Y 2n+2k+1 
Kk J 



(9) 



Dans le cas ou R(x) = x, a chaque etape de lifting, chaque 
echantillon passe-bas Y 2n peut etre exprime sous forme d'une combinaison 
lineaire d'echantillons du signal Xj, correspondant a un filtrage passe-bas du 
25 signal Xj. Ainsi, chaque etape de lifting passe-bas peut etre associee a un filtre 
passe-bas. 

On dit que chaque etape de lifting passe-bas "transforme par 
augmentation" (en anglais "lifts") un filtre passe-bas (celui qui est associe a 
I'etape de lifting passe-bas precedente) en un autre filtre passe-bas (celui qui 
30 est associe a I'etape de lifting passe-bas courante). Le terme lifting rend compte 
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du fait que le support du filtre passe-bas resultant (c'est-a-dire le nombre de 
coefficients du filtre) est plus grand que celui du filtre passe-bas correspondant 
a I'etape precedente. Le filtre passe-bas de support plus grand finalement 
obtenu par lifting est appele "filtre passe-bas equivalent". 
5 II en va de meme pour les filtres passe-haut et les etapes de lifting, 

c'est-a-dire que chaque etape de lifting passe-haut peut etre vue comme le 
lifting d'un filtre passe-haut en un filtre passe-haut de support plus grand, 
appele "filtre passe-haut equivalent". 

La meme terminologie est utilisee dans le cas ou la fonction 

10 d'approximation R n'est pas I'identite. 

Ainsi, la realisation par lifting (en anglais "lifting-based 
implementation") d'une operation de filtrage en ondelettes (H 0) Hi) peut etre vue 
comme une sequence d'operations au cours desquelles, alternativement et de 
fagon iterative, des filtres passe-bas sont transformes par augmentation en 

15 filtres passe-bas equivalents de support plus grand, et des filtres passe-haut 
sont transformes par augmentation en filtres passe-haut equivalents de support 
plus grand, ce jusqu'a ce que le filtre passe-bas soit egal a H 0 et que le filtre 
passe-haut soit egal a Hi. 

On rencontre egalement plusieurs types de symetries dans les 

20 etapes de lifting : 

- il existe deux types d'etapes de lifting symetriques : 

• des etapes de lifting a symetrie WSS et 

• des etapes de lifting a symetrie HSS ; et 

- il existe deux types d'etapes de lifting antisymetriques : 
25 • des etapes de lifting a antisymetrie WSA et 

• des etapes de lifting a antisymetrie HSA. 

Une etape de lifting symetrique de longueur impaire ou lifting WSS 
(en anglais "Whole-Sample-Symmetric lifting") est definie par une relation de 
symetrie WSS entre les coefficients de filtrage de I'etape de lifting : 
30 a_ k -a k (10) 
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Une etape de lifting symetrique de longueur paire ou lifting HSS (en 
anglais "Half-Sample-Symmetric lifting") est definie par une relation de symetrie 
HSS entre les coefficients de filtrage de I'etape de lifting : 

a_ k =a k-1 (11) 
5 Une etape de lifting antisymetrique de longueur impaire ou lifting 

WSA (en anglais "Whole-Sample-Antisymmetric lifting") est definie par une 
relation de symetrie WSA entre les coefficients de filtrage de I'etape de lifting : 

a 0 -0 et a_ k =-a k (12) 

Une etape de lifting antisymetrique de longueur paire ou lifting HSA 
10 (en anglais "Half-Sample-Antisymmetric lifting") est definie par une relation de 
symetrie HSA entre les coefficients de filtrage de I'etape de lifting : 

a_ k =-a k _-| (13) 

Une etape de lifting a coefficient unique ou lifting SC (en anglais 
"Single-Coefficient lifting") est definie par le fait que tous les coefficients de 
15 filtrage de I'etape de lifting sauf un (note a k - ) sont nuls. 

Pour les etapes de lifting passe-bas, 

- si k'=0, I'etape de lifting est appelee etape de lifting passe-bas a 
coefficient unique a droite ou lifting passe-bas RSC (en anglais "Right-side 
Single Coefficient) et 

20 - si k'= -1, I'etape de lifting est appelee etape de lifting passe-bas a 

coefficient unique a gauche ou lifting passe-bas LSC (en anglais "Left-side 
Single Coefficient). 

De fa?on analogue, pour les etapes de lifting passe-haut, 

- si k'=1, I'etape de lifting est appelee etape de lifting passe-haut a 
25 coefficient unique a droite ou lifting passe-haut RSC et 

- si k'=0, I'etape de lifting est appelee etape de lifting passe-haut a 
coefficient unique a gauche ou lifting passe-haut LSC. 

Etant donne qu'a I'heure actuelle, de plus en plus d'algorithmes de 
compression de signaux utilisent des operations de filtrage en ondelettes pour 
30 la decorrelation des signaux, il est tres interessant de disposer d'une realisation 
efficace de ces operations de filtrage. Divers algorithmes de compression 
d'images utilisent des transformations en ondelettes bidimensionnelles, qui 




consistent en I'application successive de transformations monodimensionnelles 
du type decrit en introduction. 

Une realisation efficace du filtrage monodimensionnel a base 
d'ondelettes est le filtrage a base de lifting. 
5 On sait en effet que le filtrage par lifting permet de reduire le nombre 

d'operations de filtrage (multiplications et additions) et permet egalement une 
realisation inversible du filtrage a base d'ondelettes, c'est-a-dire avec 
reconstruction parfaite du signal, sans perte d'information, dans le cas de 
signaux dont les valeurs d'echantillons sont entieres. 

10 En outre, I'etape inverse d'une etape de lifting est obtenue de fa?on 

triviale en changeant le signe devant la fonction d'approximation R donnee 
dans les equations (8) et (9) ci-dessus. Cela permet notamment d'utiliser la 
meme circuiterie dans une implantation materielle correspondant a la fois aux 
transformations en ondelettes directe (en anglais "forward 1 ) et inverse. 

15 Pour un couple de filtres en ondelettes (Ho.Hj) donne, le probleme 

consiste a trouver une sequence d'etapes de lifting qui vont realiser la 

transformation en ondelettes. C'est le probleme de !a factorisation lifting. 

On connait deux families de filtres en ondelettes ayant de bonnes 

proprietes de decorrelation de signaux : la famille WSSAA/SS, en abrege WSS, 
20 et la famille HSS/HSA, la notation A/B signifiant que le filtre passe-bas H 0 

presente une symetrie A et le filtre passe-haut Hi presente une symetrie B. 

Le probleme consistant a realiser des filtres en ondelettes WSS en 

utilisant des etapes de lifting HSS a ete resolu et est decrit dans Tarticle de 

Ingrid DAUBECHIES et Wim SWELDENS cite plus haut. 
25 Le probleme consistant a realiser des filtres en ondelettes HSS/HSA 

a ete en partie resolu, en utilisant des etapes de lifting WSA. 

Cependant, cette solution n'est que partielle puisque cette realisation 

de Tart anterieur ne s'applique qu'a un ensemble limite de filtres en ondelettes 

HSS/HSA et non a I'ensemble de ces filtres. 
30 D'autre part, la publication de E. MAJANI intitulee "Low-complexity 

wavelet filter design for image compression", TDA Progress Report 42-119, 

NASA Jet Propulsion Laboratory, Pasadena, CA, 15 novembre 1994, indique 
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qu'on peut utiliser des etapes de lifting WSA pour transformer par augmentation 
(en anglais "lift 1 ) un couple de filtres (Hq.H^), ou tes deux filtres ont le meme 

support, en un couple de filtres de support plus grand (H' 0 .H^). Cependant, on 
n'y trouve aucune indication sur la fa?on de realiser, par des etapes de lifting, 
5 des couples de filtres (Hq.H^ ou les deux filtres ont meme support et ou ce 

support est superieur a deux. 

La presente invention a pour but de remedier a cet inconvenient. 

Elle propose une realisation par lifting de tous les filtres en 
ondelettes HSS/HSA qui ne peuvent pas etre realises en utilisant uniquement 
10 des etapes de lifting WSA. 

Elle apporte done une solution au probleme de la realisation par 
lifting de couples de filtres (Hq.H^ ou les deux filtres ont le meme support et ou 

ce support est superieur a 2. 

Dans ce but, la presente invention propose un procede de filtrage 
15 adapte a transformer un signal numerique d'entree en un ou plusieurs signaux 
numeriques de sortie comportant des echantillons de rang pair et des 
echantillons de rang impair, ce procede comportant au moins une iteration qui 
contient 

- une operation de modification d'echantillons de rang pair par une 
20 fonction d'echantillons de rang impair ponderes, 

- une operation de modification d'echantillons de rang impair par une 
fonction d'echantillons de rang pair ponderes, 

les echantillons ponderes etant obtenus par au moins une operation 
de ponderation, 

25 ce procede etant remarquable en ce qu'au moins une des operations de 
ponderation est appliquee a la difference entre deux echantillons de rang pair 
consecutifs. 

Ainsi, la presente invention permet notamment la realisation par 
lifting de tout couple de filtres en ondelettes HSS/HSA, ainsi que la realisation 
30 par lifting de tout couple de filtres orthogonaux et la realisation par lifting de tout 
couple de filtres en ondelettes quelconques de support egal. 
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L'invention permet egalement, grace a la technique du filtrage par 
lifting, de reduire le nombre d'operations de filtrage (multiplications et additions) 
par rapport a un filtrage de convolution classique et garantit egalement le 
caractere inversible des transformations en ondelettes, c'est-a-dire une 
reconstruction parfaite du signal, dans le cas de signaux dont les valeurs 
d'echantillons sont entieres. 

Selon une caracteristique particuliere, 1'operation de modification 
d'echantillons de rang impair est effectuee a la suite de 1'operation de 
modification d'echantillons de rang pair. 

Dans un mode particulier de realisation, I'iteration precitee consiste 
notamment a : 

- ponderer, par un premier coefficient de ponderation, au moins un 
echantillon de rang impair voisin d'un echantillon pair en cours de modification, 
de fagon a obtenir un echantillon de rang impair pondere, 

- modifier au moins un echantillon de rang pair a partir d'au moins 
un echantillon de rang impair pondere, 

- ponderer, par un deuxieme coefficient de ponderation, des. 
echantillons de rang pair voisins d'un echantillon impair en cours de 
modification, de fagon a obtenir des echantillons de rang pair ponderes, et 

- modifier au moins un echantillon de rang impair a partir d'au moins 
un echantillon de rang pair pondere. 

Par echantillon "voisin", on entend un echantillon dont le rang est 
consecutif a I'echantillon considere. 

Selon une caracteristique particuliere, le deuxieme coefficient de 
ponderation est fonction du premier coefficient de ponderation. 

Le deuxieme coefficient de ponderation peut par exemple dependre 
du premier coefficient de ponderation comme suit : 



Poj =n V 



1-2la 0ii 
i=0 



ou <xo,i designe le premier coefficient de ponderation, p 0 ,j designe le deuxieme 
coefficient de ponderation, i et j sont des entiers et m, est une valeur definie par 

la recurrence m 0 = (-1) L ° et nrij = -m H , L 0 etant un entier predetermine. 

Selon une caracteristique particuliere, a chaque iteration, 
5 Techantillon de rang impair voisin de Techantillon pair en cours de modification 
est alternativement Techantillon de rang immediatement inferieur ou 
immediatement superieur. 

Selon une caracteristique particuliere, le procede de filtrage 
comporte, a Tissue de Titeration precitee, une etape de filtrage supplemental 
10 incluant une operation de ponderation par un troisieme coefficient de 
ponderation. 

Selon une caracteristique particuliere, le troisieme coefficient de 
ponderation est fonction du coefficient de ponderation utilise a Tetape 
precedente, comme suit : 
15 Y = -1'(2P 0l L 0 -l) 

ou y designe le troisieme coefficient de ponderation, L 0 est un parametre 
predetermine et Po,l 0 -1 d6s '9 ne le coefficient de ponderation utilise a Tetape 
precedente. 

Dans ('application de ('invention au filtrage d'images numeriques, le 
20 signal numerique d'entree represente une image. 

Selon un deuxieme aspect, la presente invention propose aussi un 
procede de filtrage adapte a transformer un ou plusieurs signaux numeriques 
d'entree en un signal numerique de sortie, les signaux d'entree comportant des 
echantillons de rang pair et des echantillons de rang impair, le procede 
25 comportant au moins une iteration qui contient 

- une operation de modification d'echantillons de rang impair par une 
fonction d'echantillons de rang pair ponderes, 

- une operation de modification d'echantillons de rang pair par une 
fonction d'echantillons de rang impair ponderes, 

30 les echantillons ponderes etant obtenus par au moins une operation 

de ponderation, 
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ce procede etant remarquable en ce qu'au moins une des operations de 
ponderation est appliquee a la difference entre deux echantillons de rang pair 
consecutifs. 

Selon une caracteristique particuliere, I'operation de modification 
5 d'echantillons de rang pair est effectuee a la suite de I'operation de modification 
d'echantillons de rang impair. 

Dans un mode particulier de realisation, I'iteration precitee consiste 
notamment a : 

- ponderer, par un quatrieme coefficient de ponderation, des 
10 echantillons de rang pair voisins d'un echantillon impair en cours de 

modification, de fa?on a obtenir des echantillons de rang pair ponderes, 

- modifier au moins un echantillon de rang impair a partir d'au moins 
un echantillon de rang pair pondere, 

- ponderer, par un cinquieme coefficient de ponderation, au moins 
15 un echantillon de rang impair voisin d'un echantillon pair en cours de 

modification, de fa?on a obtenir un echantillon de rang impair pondere et 

- modifier au moins un echantillon de rang pair a partir d'au moins 
un echantillon de rang impair pondere. 

Selon une caracteristique particuliere, le quatrieme coefficient de 
20 ponderation est fonction du cinquieme coefficient de ponderation. 

Le quatrieme coefficient de ponderation peut par exemple dependre 
du cinquieme coefficient de ponderation comme suit : 



' J ^ 
1-25> 0i , 
i=0 



ou oco.i designe le cinquieme coefficient de ponderation, p 0 j designe le 
25 quatrieme coefficient de ponderation, i et j sont des entiers et mj est une valeur 

definie par la recurrence m 0 =(-1) L ° et mj = -rrij-i, L 0 etant un entier 

predetermine. 

Suivant le deuxieme aspect de ('invention, selon une caracteristique 
particuliere, a chaque iteration, I'echantillon de rang impair voisin de 
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I'echantillon pair en cours de modification est alternativement I'echantillon de 
rang immediatement inferieur ou immediatement superieur. 

Suivant le deuxieme aspect de I'invention, selon une caracteristique 
particuliere, le procede de filtrage comporte, a Tissue de Iteration precitee, une 
5 etape de filtrage supplementaire incluant une operation de ponderation par un 
sixieme coefficient de ponderation. 

Selon une caracteristique particuliere, le sixieme coefficient de 
ponderation est fonction du coefficient de ponderation utilise a I'etape 
precedente, comme suit : 

10 Y = -1/t2p 0 ,L 0 -l) 

ou y designe le sixieme coefficient de ponderation, L 0 est un parametre 
predetermine et Po,l 0 -i designe le coefficient de ponderation utilise a I'etape 

precedente. 

Dans I'application de I'invention au filtrage d'images numeriques, le 
15 signal numerique de sortie represente une image. 

Selon une caracteristique particuliere des premier et second aspects 
de I'invention, les operations de modification consistent a appliquer une fonction 
d'approximation. 

Cette fonction peut etre : 
20 - la fonction identite, ou 

- une fonction d'une variable reelle qui fournit rentier le plus proche 
de la variable, ou 

- une fonction d'une variable reelle qui fournit le premier entier 
inferieur a la variable, ou 

25 - une fonction d'une variable reelle qui fournit le premier entier 

superieur a la variable, ou encore 

- une fonction d'une variable decomposee en sous-variables dont la 
somme est egale a la variable, qui fournit une somme de valeurs approchees 
des sous-variables, chacune des valeurs approchees des sous-variables etant, 

30 soit une fonction d'une variable reelle qui fournit rentier le plus proche de la 
variable, soit une fonction d'une variable reelle qui fournit le premier entier 
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inferieur a la variable, soit une fonction d'une variable reelle qui fournit le 
premier entier superieur a la variable. 

L'invention vise aussi un dispositif de traitement de signal qui 
comporte des moyens adaptes a mettre en oeuvre un precede de filtrage tel que 
ci-dessus. 

Dans le meme but que celui indique plus haut, la presente invention 
propose egalement, selon un troisieme aspect, un dispositif de filtrage 
numerique adapte a transformer un signal numerique d'entree en un ou 
plusieurs signaux numeriques de sortie comportant des echantillons de rang 
pair et des echantillons de rang impair, ce filtre comportant 

- au moins un module de ponderation, 

- un module pour modifier des echantillons de rang pair par une 
fonction d'echantillons de rang impair ponderes, 

- un module pour modifier des echantillons de rang impair par une 
fonction d'echantillons de rang pair ponderes, 

ces echantillons ponderes etant fournis par le module de 
ponderation, les modules de modification fonctionnant de fa?on iterative, de 
fa?on a modifier au moins une fois des echantillons de rang pair puis au moins 
une fois des echantillons de rang impair, 

ce dispositif de filtrage etant remarquable en ce qu'au moins un des modules de 
ponderation precites re9oit en entree la difference entre deux echantillons de 
rang pair consecutifs. 

Dans le meme but que celui indique plus haut, la presente invention 
propose egalement, selon un quatrieme aspect, un dispositif de filtrage 
numerique adapte a transformer un ou plusieurs signaux numeriques d'entree 
en un signal numerique de sortie, les signaux d'entree comportant des 
echantillons de rang pair et des echantillons de rang impair, ce dispositif de 
filtrage comportant 

- au moins un module de ponderation, 

- un module pour modifier des echantillons de rang impair par une 
fonction d'echantillons de rang pair ponderes, 
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- un module pour modifier des echantillons de rang pair par une 
fonction d'echantillons de rang impair ponderes, 

les echantillons ponderes etant fournis par le module de ponderation, 
le module de modification fonctionnant de fa?on iterative, de fa<?on a modifier au 
5 moins une fois des echantillons de rang impair puis au moins une fois des 
echantillons de rang pair, 

ce dispositif de filtrage etant remarquable en ce qu'au moins un des modules de 
ponderation precites re?oit en entree la difference entre deux echantillons de 
rang pair consecutifs. 

10 Les caracteristiques particulieres et les avantages du dispositif de 

filtrage etant similaires a ceux du procede de filtrage selon I'invention, ils ne 
sont pas rappeles ici. 

Toujours dans le meme but, I'invention propose aussi un dispositif de 
traitement de signal incluant un dispositif de filtrage tel que precedemment 
15 defini, ou les moyens de mise en oeuvre du procede correspondant. 

L'invention concerne aussi un dispositif de traitement de signal 
incluant au moins deux dispositifs de filtrage tels que precedemment definis, le 
signal de sortie de Tun des dispositifs de filtrage etant le signal d'entree de 
I'autre dispositif de filtrage. 
20 L'invention concerne aussi un appareil numerique incluant le 

dispositif de traitement de signal. 

L'invention concerne aussi un appareil photographique numerique 
incluant le dispositif de traitement de signal. 

Toujours dans le meme but, la presente invention propose aussi un 
25 procede de codage comportant des etapes adaptees a mettre en oeuvre un 
procede de filtrage tel que ci-dessus. 

Toujours dans le meme but, la presente invention propose aussi un 
dispositif de codage comportant au moins un dispositif de filtrage tel que ci- 
dessus, ou les moyens de mise en oeuvre du procede correspondant. 
30 Toujours dans le meme but, la presente invention propose aussi un 

procede de compression de signaux numeriques, comportant des etapes 
adaptees a mettre en oeuvre un procede de filtrage tel que ci-dessus. 
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Toujours dans le meme but, la presente invention propose aussi un 
dispositif de compression de signaux numeriques, comportant au moins un 
dispositif de filtrage tel que ci-dessus, ou les moyens de mise en oeuvre du 
procede correspondant. 
5 Un moyen de stockage d'information, lisible par un ordinateur ou par 

un microprocesseur, integre ou non au dispositif, eventuellement amovible, 
memorise un programme mettant en oeuvre un procede de filtrage tel que ci- 
dessus. 

La presente invention vise aussi un produit programme d'ordinateur 
10 comportant des sequences d'instructions pour mettre en oeuvre un procede de 
filtrage tel que ci-dessus. 

Les caracteristiques particulieres et les avantages des dispositifs de 
traitement de signal, de I'appareil numerique, de I'appareil photographique 
numerique, des procedes et dispositifs de codage et de compression, du moyen 
15 de stockage et du produit programme d'ordinateur etant similaires a ceux du 
procede de filtrage selon I'invention, ils ne sont pas rappeles ici. 

D'autres aspects et avantages de ('invention apparaTtront a la lecture 
de la description detaillee qui suit d'un mode particulier de realisation, donne a 
titre d'exemple non limitatif. La description se refere aux dessins qui 
20 I'accompagnent, dans lesquels : 

- la figure 1 represente de maniere schematique un dispositif de 
traitement de donnees selon Invention, 

- la figure 2 represente de maniere schematique un autre dispositif 
de traitement de donnees selon I'invention, 

25 - la figure 3 represente un mode de realisation du dispositif de 

traitement de donnees selon I'invention, 

- la figure 4 represente un mode de realisation de circuit de 
transformation inclus dans le dispositif de traitement de donnees de la figure 1, 

- la figure 5 est un organigramme illustrant les principales etapes 
30 d'un procede de filtrage conforme a la presente invention, dans un premier 

mode particulier de realisation ; 




- la figure 6 est un organigramme illustrant les principales etapes 
d'un procede de filtrage conforme a la presente invention, dans un deuxieme 
mode particulier de realisation ; 

- la figure 7 est un organigramme illustrant un mode de realisation 
5 de la determination des parametres de lifting a partir des coefficients des filtres 

en ondelettes ; 

- les figures 8A, 8B et 8C sont des organigrammes illustrant 
respectivement les procedures de calcul notees UNLIFT_H0, UNLIFT_H1 et 
UNLIFT_EQ intervenant dans I'organigramme de la figure 7 ; 

10 - la figure 9 represente un premier mode de realisation d'un bloc de 

transformation elementaire selon la presente invention ; et 

- la figure 10 represente un deuxieme mode de realisation d'un bloc 
de transformation elementaire selon la presente invention. 

Selon un mode de realisation choisi et represente a la figure 1, un 

15 dispositif de traitement de donnees selon I'invention est un dispositif 2 de 
codage de donnees qui comporte une entree 24 a laquelle est reliee une source 
1 de donnees non codees. Le dispositif de traitement de donnees peut etre 
integre a un appareil numerique, tel qu'un ordinateur, un appareil photo 
numerique, un telecopieur, parexemple. 

20 La source 1 comporte par exemple un moyen de memoire, tel que 

memoire vive, disque dur, disquette, disque compact, pour memoriser des 
donnees non codees, ce moyen de memoire etant associe a un moyen de 
lecture approprie pour y lire les donnees. Un moyen pour enregistrer les 
donnees dans le moyen de memoire peut egalement etre prevu. La source 1 

25 peut egalement etre integree a I'appareil numerique. 

On considerera plus particulierement dans la suite que les donnees a 
coder sont une suite d'echantillons numeriques representant une image IM. 

La source 1 fournit un signal numerique d'image SI a I'entree du 
circuit de codage 2. Le signal d'image SI est une suite de mots numeriques, par 

30 exemple des octets. Chaque valeur d'octet represente un pixel de I'image IM, ici 
a 256 niveaux de gris, ou image noir et blanc. L'image peut etre une image 
multispectrale, par exemple une image en couleur ayant des composantes dans 
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trois bandes de frequence, de type rouge-vert-bleu ou luminance et 
chrominance. Chaque bande est alors traitee de maniere analogue a I'image 
monospectrale. 

Des moyens 3 utilisateurs de donnees codees sont relies en sortie 
5 25 du dispositif de codage 2. 

Les moyens utilisateurs 3 comportent par exemple des moyens de 
memorisation de donnees codees, et/ou des moyens de transmission des 
donnees codees. 

Le dispositif de codage 2 comporte classiquement, a partir de 
10 I'entree 24, un circuit de transformation 21, plus particulierement concerne par 
la presente invention, et dont un exemple de realisation sera detaille dans la 
suite. Les transformations envisagees ici sont des decompositions en sous- 
bandes de frequence du signal de donnees, de maniere a effectuer une analyse 
du signal. 

15 Le circuit de transformation 21 est relie a un circuit de quantification 

22. Le circuit de quantification met en oeuvre une quantification connue en soi, 
par exemple une quantification scalaire, ou une quantification vectorielle, des 
coefficients, ou de groupes de coefficients, des sous-bandes de frequence 
fournies par le circuit 21 . 

20 Le circuit 22 est relie a un circuit 23 de codage entropique, qui 

effectue un codage entropique, par exemple un codage de Huffman, ou un 
codage arithmetique, des donnees quantifies par le circuit 22. 

La figure 2 represente un autre dispositif de traitement de donnees 
selon I'invention, sous la forme d'un dispositif 5 de decodage de donnees 

25 codees par le dispositif 2. 

Des moyens 4 utilisateurs de donnees codees sont relies en entree 

54 du dispositif de decodage 5. Les moyens 4 comportent par exemple des 
moyens de memoire de donnees codees, et/ou des moyens de reception des 
donnees codees qui sont adaptes a recevoir les donnees codees transmises 

30 par les moyens de transmission 3. 

Des moyens 6 utilisateurs de donnees decodees sont relies en sortie 

55 du dispositif de decodage 5. Les moyens utilisateurs 6 sont par exemple des 
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moyens de visualisation d'images, ou des moyens de restitution de sons, en 
fonction de la nature des donnees traitees. 

Le dispositif de decodage 5 effectue globalement des operations 
inverses de celles du dispositif de codage 2. Le dispositif 5 comporte un circuit 
5 de decodage entropique 51, qui effectue un decodage entropique 
correspondant au codage du circuit 23. Le circuit 51 est relie a un circuit de 
dequantification 52, correspondant au circuit de quantification 22. 

Le circuit 52 est relie a un circuit de transformation inverse 53, 
correspondant au circuit de transformation 21. Le circuit de transformation 
10 inverse 53 est plus particulierement concerne par la presente invention. Un 
exemple de realisation sera detaille dans la suite. Les transformations 
envisagees ici realisent une synthese du signal numerique, a partir de sous- 
bandes de frequence. 

Le dispositif de codage et le dispositif de decodage peuvent etre 
15 integres dans un meme appareil numerique, par exemple une camera 
numerique. Dans ce cas, le dispositif de traitement de donnees realise le 
codage et le decodage des donnees. 

En reference a la figure 3, est decrit un exemple de dispositif 10 
mettant en oeuvre invention. Ce dispositif est adapte a transformer un signal 
20 numerique, et selon les exemples developpes dans la suite, a I'analyser, ou a le 
synthetiser ou encore a I'analyser puis le synthetiser. Ainsi, la figure 3 illustre un 
dispositif de codage ou de decodage. 

Le dispositif 10 est ici un micro-ordinateur comportant un bus de 
communication 101 auquel sont relies : 
25 - une unite centrale de traitement 100, 

- une memoire morte 102, 

- une memoire vive 103, 

- un ecran 104, 

- un clavier 114, 

30 - un disque dur 108, 

- un lecteur de disquettes 109 adapte a recevoir urje disquette 110, 
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- une interface 112 de communication avec un reseau de 
communication 113, 

- une carte d'entrees/sorties 106 reliee a un microphone 111. 

La memoire vive 103 conserve des donnees, des variables et des 
5 resultats intermediates de traitement, dans des registres de memoire portant, 
dans la description, les memes noms que les donnees dont ils conservent les 
valeurs. La memoire vive 103 comporte notamment : 

- un registre "/c", contenant la valeur courante de la variable k definie 

plus loin, 

10 - un registre "/', contenant la valeur courante de la variable j definie 

plus loin, 

- un registre "y/\ dans lequel sont conservees les valeurs des 
echantillons yj de I'image transformee au fur et a mesure de leur calcul, 

- un registre "fikj\ dans lequel sont conservees les valeurs des 
15 coefficients de ponderation p k j, definis plus loin, 

- un registre "m", dans lequel sont conservees les valeurs des 
parametres rrij, definis plus loin, 

- un registre "f, dans lequel sont conservees les valeurs des 
coefficients y, definis plus loin. 

20 Dans le cas ou la figure 3 illustre un dispositif de codage, la memoire 

morte 102 est adaptee a conserver notamment : 

- le programme de fonctionnement de I'unite centrale de traitement 
100, dans un registre "Program", 

- la valeur maximale k max de la variable k, dans un registre "k max % \ 

25 - la valeur des parametres L kl definis plus loin, pour toutes les 

valeurs de la variable k comprises entre 0 et k max , dans un registre "L k ", 

- les valeurs des coefficients de ponderation a k ,j, dans un registre 

a Ki • 

II est a noter que, dans le cas ou la figure 3 illustre un dispositif de 
30 decodage, le registre "Program" est egalement en memoire morte mais les 
registres "k max " % "L k " et "a kJ " sont en memoire vive. 
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L'unite centrale de traitement 100 est adaptee a mettre en oeuvre les 
organigrammes des figures 5 et 7. 

Le disque dur 108 memorise les programmes mettant en oeuvre 
Tinvention, ainsi que les donnees a coder et les donnees codees selon 
5 Tinvention. Ces programmes peuvent aussi etre lus sur la disquette 110, ou 
regus via le reseau de communication 113, ou encore memorises en memoire 
morte 102. 

De maniere plus generale, les programmes selon la presente 
invention sont memorises dans un moyen de stockage. Ce moyen de stockage 

10 est lisible par un ordinateur ou par un microprocesseur. Ce moyen de stockage 
est integre ou non au dispositif, et peut etre amovible. Par exemple, il peut 
comporter une bande magnetique, une disquette ou un CD-ROM (disque 
compact a memoire figee). 

Le dispositif 10 peut recevoir des donnees a coder en provenance 

15 d'un dispositif peripherique 107, tel qu'un appareil photographique numerique, 
ou un scanner, ou tout autre moyen d'acquisition ou de stockage de donnees. 

Le dispositif 10 peut egalement recevoir des donnees a coder en 
provenance d'un dispositif distant, via le reseau de communication 113, et 
transmettre des donnees codees vers un dispositif distant, toujours via le 

20 reseau de communication 113. 

Le dispositif 10 peut aussi recevoir des donnees a coder en 
provenance du microphone 111. Ces donnees sont alors un signal de son. 

L'ecran 104 permet a un utilisateur notamment de visualiser les 
donnees a coder, et sert, avec le clavier 1 14, d'interface utilisateur. 

25 En reference a la figure 4, ie circuit de transformation 21, ou circuit 

d'analyse, est un circuit de decomposition dyadique a deux niveaux. Le circuit 
21 est, dans ce mode de realisation, un ensemble classique de filtres, 
respectivement associes a des decimateurs par deux, qui filtrent le signal 
d'image selon deux directions respectivement verticale et horizontale, en 

30 signaux de sous-bandes de hautes et basses frequences spatiales. La relation 
entre un filtre passe-haut et un filtre passe-bas est determinee par les 
conditions de reconstruction parfaite du signal. Differents exemples de filtres 
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vont etre envisages dans la suite. II est a noter que les filtres de decomposition 
verticale et horizontale ne sont pas forcement identiques, bien qu'en pratique ce 
soit generalement le cas. Le circuit 21 comporte ici deux blocs successifs 
d'analyse B1 et B2 pour decomposer une image IM en des signaux de sous- 
5 bandes selon deux niveaux de resolution. 

De maniere generate, la resolution d'un signal est le nombre 
d'echantillons par unite de longueur utilises pour representer ce signal. Dans le 
cas d'un signal d'image, la resolution d'un signal de sous-bande est liee au 
nombre d'echantillons par unite de longueur utilises pour representer ce signal 

10 de sous-bande horizontalement et verticalement. La resolution depend du 
nombre de decimations effectuees, du facteur de decimation et de la resolution 
de I'image initiale. 

Le premier bloc d'analyse B1 retpoit un signal numerique d'image et 
I'applique a deux filtres numeriques respectivement passe-bas et passe-haut 

15 210 et 220 qui filtrent le signal d'image selon une premiere direction, par 
exemple horizontale dans le cas d'un signal d'image. Apres passage par des 
decimateurs par deux D210 et D220, les signaux filtres resultants sont 
respectivement appliques a deux filtres passe-bas 230 et 250, et passe-haut 
240 et 260, qui les filtrent selon une seconde direction, par exemple verticale 

20 dans le cas d'un signal d'image. Chaque signal filtre resultant passe par un 
decimateur par deux respectif D230, D240, D250 et D260. Le bloc B1 delivre en 
sortie quatre signaux de sous-bandes LL^ LHi, HLi et Hhh de resolution RESi 
la plus elevee dans la decomposition. 

Le signal de sous-bande LLi comporte les composantes, ou 

25 coefficients, de basse frequence, selon les deux directions, du signal d'image. 
Le signal de sous-bande LHi comporte les composantes de basse frequence 
selon une premiere direction et de haute frequence selon une seconde 
direction, du signal d'image. Le signal de sous-bande HLt comporte les 
composantes de haute frequence selon la premiere direction et les 

30 composantes de basse frequence selon la seconde direction. Enfin, le signal de 
sous-bande Hhh comporte les composantes de haute frequence selon les deux 
directions. 
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Chaque signal de sous-bande est un ensemble de coefficients reels 
constant a partir de I'image d'origine, qui contient de reformation correspondant 
a une orientation respectivement verticale, horizontale et diagonale des 
contours de Timage, dans une bande de frequence donnee. Chaque signal de 
5 sous-bande peut etre assimile a une image. 

Le signal de sous-bande est analyse par le bloc d'analyse B2 
analogue au precedent pour fournir quatre signaux de sous-bandes LL 2 , LH 2 , 
HL 2 et HH 2 de niveau de resolution RES 2 . 

Chacun des signaux de sous-bandes de resolution RES 2 correspond 
10 egalement a une orientation dans I'image. 

A un circuit d'analyse 21 donne, correspond, de maniere classique, 
un circuit de synthese dont la structure se deduit de celle du circuit d'analyse. 

La figure 5 montre la succession des etapes de lifting realisant la 
transformation en ondelettes directe conformement a la presente invention. La 
15 transformation inverse sera decrite ensuite en liaison avec la figure 6. 

On considere un signal d'entree x n a transformer, ayant des 
echantillons dont les rangs sont compris entre deux entiers i 0 et h (ii non 
compris), en un signal de sortie y n ayant des echantillons dont les rangs sont 
compris entre i 0 et u (ii non compris). 
20 Une etape d'initialisation 500 consiste tout d'abord a mettre la valeur 

des echantillons du signal de sortie y n a celle des echantillons du signal 
d'entree x n . 

Lors de I'etape 500, on initialise egalement deux compteurs, dont la 
valeur courante sera notee respectivement k et j, a la valeur 0. 
25 Lors de I'etape 500, on initialise aussi le parametre m,, mentionne 

plus haut en relation avec la figure 3,. a une valeur m 0 =(-1) L ° , ou L 0 est un 

autre parametre, defini plus loin. 

Si on considere un signal y n ayant des echantillons dont les rangs 
sont compris entre i 0 et ii (ii non compris), le signal y n doit etre etendu au-dela 
30 de cette plage, suffisamment pour que tous les echantillons du signal de sortie 
y n dont les rangs sont compris entre i 0 et h (h non compris) puissent etre 
calcules. Cela ressortira de la description d'un exemple, donnee plus loin. 
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Ainsi, une procedure d'extension de signal (non representee sur la 
figure 5) consiste a calculer les valeurs du signal de sortie yj pour des valeurs 
de i sortant d'un intervalle predetermine i 0 < i < i^ : 

Yi =ypSE(i.i 0 .ii) ( 14 ) 
5 ou PSE(Uo.ii)) = io + min[mod(i-i 0) 2(i 1 -io)) ) 2(ii-io)-mod(i-io,2(i 1 -i 0 ))] (1 5) 

ou min(a.b) d^signe la fonction fournissant la plus petite des valeurs a et b et 
mod(a.b) designe rentier compris entre 0 et b-1 (b-1 inclus) egal a la valeur de 
a, a un multiple de b pres. 

L'exemple, nullement limitatif, de procedure d'extension donne ci- 
10 dessus s'applique a des signaux tels que les rangs i 0 et h sont pairs. D'autres 
procedures d'extension peuvent etre utilisees dans le cas ou Tun au moins des 
rangs io et ii est impair. Neanmoins, le reste de la description s'applique a tous 
les cas de parite. Par ailleurs, il existe d'autres procedures d'extension pour 
d'autres types de filtres, comme par exemple, I'extension dite circulaire, connue 
15 de I'homme du metier, pour les filtres orthogonaux. 

L'etape suivante 502 du procede de filtrage est la premiere etape de 
lifting. C'est une etape de lifting passe-haut. 

Dans un mode de realisation prefere, cette etape de lifting est du 

type, soit LSC, soit RSC, en fonction de la valeur du parametre m 0 = (-1) L ° , et 
20 est definie par : 

Vn,y 2n +1 = Y2n + 1 -R(y2n+1-m 0 ) ( 16 ) 
Cette etape de lifting correspond au filtre passe-haut HSA 
H-|(z) = z-1 ou 1-z, en fonction de la valeur de m 0 . Elle correspond a une 
augmentation d'une unite du support du filtre passe-haut equivalent. 
25 A Tissue de l'etape de lifting 502, on effectue un test 504 pour 

determiner si la valeur du compteur j est egale ou non a L 0 . 

Tant que le resultat du test 504 est negatif, on itere une etape 506 
qui consiste essentiellement en une sequence d'operations de modification et 
de ponderation des echantillons de rang pair ou impair du signal y n . 
30 En effet, lors de l'etape 506, on effectue une succession de L 0 

couples d'etapes de lifting, correspondant a des valeurs du compteur j 
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comprises entre 0 et L 0 -1, ou la premiere etape de lifting est une etape de lifting 
passe-bas a coefficient unique : 

Vn,y 2 n =y2n+R(aoj-y2n+mj) (17) 
et ou la seconde etape de lifting est une etape de lifting passe-haut HSA : 
5 Vn,y 2n+ 1 =y2n+1+R(Po l j-(y2n -y2n+2» ( 18 > 

Le parametre non nul a 0 j utilise dans I'equation (17) est une 
fonction des coefficients du couple de filtres en ondelettes (Hq.Hi) et son calcul 

est detaille plus loin. 

Etant donne que la valeur de rrij est soit -1, soit +1 eomme on va le 
10 voir, le rang de I'echantillon y2n+mj est so 't 2n-1, soit 2n+1, c'est-a-dire un 

entier consecutif au rang de I'echantillon y 2n en cours de modification. On dit 
que I'echantillon de rang impair y2n+mj est voisin de I'echantillon y2n- 

L'etape de lifting passe-bas a coefficient unique est du type, soit 
LSC, soit RSC, en fonction de la valeur du parametre mj, qui est defini de fagon 
15 recursive par : rrij = -mj_ 1 . 

II y a alternance systematique des etapes suivantes (c'est-a-dire 
correspondant a des valeurs suivantes de j) de lifting passe-bas a coefficient 
unique : par exemple, la deuxieme est du type LSC si la premiere est du type 
RSC ou vice versa. La derniere etape de lifting passe-bas a coefficient unique 

20 est toujours du type LSC. 

L'etape de lifting passe-bas de ('equation (17) garantit que le support 
du filtre passe-bas equivalent augmente jusqu'a atteindre la valeur du support 
du filtre passe-haut equivalent de l'etape de lifting passe-haut precedente, i.e. la 
valeur paire superieure la plus proche. 

25 Conformement a un mode de realisation prefere de la presente 

invention, l'etape de lifting passe-haut de I'equation (18) est du type HSA, 
contrairement a I'art anterieur, qui utilise exclusivement des etapes de lifting a 
coefficient unique et des etapes de lifting WSA pour la realisation par lifting de 
filtres en ondelettes HSS/HSA. 
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L'etape de lifting passe-haut de ('equation (18) garantit que le support 
du filtre passe-haut equivalent augmente d'une valeur de deux (c'est-a-dire 
jusqu'a atteindre la valeur paire superieure la plus proche) par rapport au 
support du filtre passe-haut equivalent de l'etape de lifting passe-haut 
5 precedente. 

La filtre passe-haut equivalent de l'etape de lifting passe-haut de 
I'equation (18) est du type HSA si et seulement si le filtre passe-haut equivalent 
de l'etape de lifting passe-haut precedente est HSA et le parametre non nul Pq j 

utilise dans I'equation (18) est une fonction du parametre mj et des parametres 

10 ao ( j provenant des etapes de lifting passe-bas a coefficient unique 

precedentes, cette fonction etant par exemple definie comme suit : 

( J ^ 



Poj =*V 



(19) 



A la difference de l'etape de lifting passe-bas de I'equation (17), qui 
transforme par augmentation les coefficients passe-bas a partir des coefficients 
15 passe-haut alternativement a gauche et a droite, l'etape de lifting passe-haut de 
I'equation (18) transforme par augmentation les coefficients passe-haut a partir 
des coefficients passe-bas a la fois a gauche et a droite. 

Lorsque le resultat du test 504 est positif, c'est-a-dire que j = L 0 , 
I'alternance d'etapes de lifting passe-bas SC et passe-haut HSA prend fin et est 
20 suivie d'une etape de lifting passe-bas RSC 508, definie par : 

vn.y 2 n=y2n+R(y-y2n + i) (20) 

Dans un mode de realisation prefere, le coefficient y est calcule en 
fonction du coefficient de ponderation utilise a l'etape de lifting passe-haut 
precedente : 

25 Y = -1'( 2 Po,L 0 -l) < 21 > 

afin de garantir que le filtre passe-bas correspondant soit HSS. 

L'etape de lifting passe-bas 508 garantit que le support du filtre 
passe-bas equivalent augmente jusqu'a atteindre la valeur du support du filtre 
passe-haut equivalent de l'etape de lifting passe-haut precedente, c'est-a-dire, 



25 



dans le mode prefere de realisation decrit ici, rentier pair superieur le plus 
proche. 

Dans ce mode prefere de realisation, on obtient, a Tissue de I'etape 
508, un couple correspondant de filtres passe-bas et passe-haut qui sont de 
5 support pair egal et qui sont respectivement HSS et HSA. 

On procede ensuite a une iteration d'etapes de lifting passe-bas 
WSA (etape 512) et passe-haut WSA (etape 518), soit k max -1 etapes de 

lifting WSA au total, la valeur du compteur k croissant de 1 a k max -1. 

Pour cela, a la suite de I'etape 508, on incremente d'une unite la 
10 valeur d'un compteur k a I'etape 510. On incremente a nouveau la valeur de k 
apres chaque passage par I'etape 512 (voir etape 514 sur la figure 5) et apres 
chaque passage par I'etape 518 (voir etape 520 sur la figure 5). 

On reitere les etapes 512, 514, 518 et 520 jusqu'a ce que la valeur 
de k soit egale a k max (voir tests 516 et 522 sur la valeur de k, respectivement 
15 effectues apres les etapes decrementation 514 et 520). 

Les etapes de lifting passe-bas WSA 512 sont definies par : 

f«-2k-1-1 ^ 



Vn,y 2n =y2n + R 



X a 2k-1,j-(y2n + 1-2j -y2n+1+2j) 

ho 



(22). 



et les etapes de lifting passe-haut WSA 518 sont definies par : 



Vn,y 2n + 1 = Y2n + 1 +R 



S a 2k,j (y2n-2j ~y2n + 2j) 
j = 0 j 



(23) 



20 Ces etapes sont connues en soi. Elles ont pour effet une 

augmentation de la taille du filtre correspondant (passe-bas pour I'etape 512 et 
passe-haut pour I'etape 518) d'un multiple de quatre. 

Pour realiser des filtres en ondelettes avec un gain en courant 
continu ou gain DC (en anglais "Direct Current') unite et un gain de Nyquist de 

25 2, c f est-a-dire avec une normalisation dite (1,2), aucune etape specifique de 
normalisation du type yj = G.yj, ou G est une valeur reelle, n'est necessaire. En 
revanche, si une normalisation differente est requise, telle qu'une normalisation 
(Gi, 2G 2 ) avec un gain DC de et un gain de Nyquist de 2G 2l une etape finale 
de normalisation est necessaire : y 2n =Gj.y 2 n et y2n+i = G 2 y2n+i- 
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La fonction d'approximation R n'est pas necessairement la meme 
pour toutes les etapes de lifting dans lesquelles elle est utilisee. 

Les parametres a choisir pour la realisation lifting des filtres 
HSS/HSA sont L k ,a k j, k max et la fonction R a chaque etape de lifting. On peut 

ainsi etablir une mise en correspondance (en anglais "one-to-one mapping") 
entre ces parametres et les coefficients et longueurs de tout couple de filtres en 
ondelettes HSS/HSA a reconstruction parfaite. 

La figure 6 montre la succession des etapes de lifting realisant la 
transformation en ondelettes inverse conformement a la presente invention. 

On considere un signal d'entree y n a transformer, ayant des 
echantillons dont les rangs sont compris entre deux entiers i 0 et ii (h non 
compris) en un signal de sortie x n ayant des echantillons dont les rangs sont 
compris entre io et ii (ii non compris). 

Une etape d'initialisation 600 consiste tout d'abord a mettre la valeur 
des echantillons du signal de sortie x n a celle des echantillons du signal 
d'entree y n . 

Toutes les etapes de la realisation lifting inverse sont les operations 
inverses des etapes de lifting effectuees dans la transformation directe. Par 
ailleurs, les etapes de la realisation lifting inverse sont effectuees dans I'ordre 
contraire des etapes de la realisation lifting directe. 

Pour un traitement d'analyse suivi d'un traitement de synthese des 
memes signaux, tous les parametres de la transformation inverse sont 
identiques a ceux choisis pour la transformation directe. 

Comme le montre la figure 6, lors de I'etape d'initialisation 600, on 
initialise egalement un compteur, dont la valeur courante est notee k, a la valeur 

^max ~ 1 • 

Si on considere un signal x n ayant des echantillons dont les rangs 
sont compris entre i 0 et h (\^ non compris), le signal x n doit etre etendu au-dela 
de cette plage, suffisamment pour que tous les echantillons du signal de sortie 
x n dont les rangs sont compris entre i 0 et ii (h non compris) puissent etre 
calcules. Cela ressortira de la description d'un exemple, donnee plus loin. 
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Ainsi, une procedure d'extension de signal (non representee sur la 
figure 6) consiste a calculer les valeurs du signal de sortie x s pour des valeurs 
de i sortant d'un intervalle predetermine \q < i < if, 

- comme defini par I'equation (24) ci-dessous pour des valeurs de i 
5 telles que mod(i-i 0l 2(ii-io)) < Mo : 

x i = x i 0 +mod(i-i 0t 2(i 1 -i 0 )) ( 24 ) 
ou mod(a.b) designe rentier compris entre 0 et b-1 (b-1 inclus) egal a la valeur 
de a, a un multiple de b pres, 

- comme defini par I'equation (25) ci-dessous pour des valeurs 
10 paires de i telles que mod(i-i 0) 2(ii-i 0 )) > Mo : 

x i = x i 0 + mod(2(i 1 -i 0 -1)-(i-i 0 ),2(ii -i 0 ))) ( 25 ) 

- et comme defini par I'equation (26) ci-dessous pour des valeurs 
impaires de i telles que mod(i-i 0l 2(Mo)) > Mo : 

x i = x i 0 +mod(2(i 1 - i 0 )-(i-i 0 ).2(h -i 0 ))) < 26 > 
15 II existe d'autres procedures d'extension pour d'autres types de 

filtres, comme par exemple, I'extension dite circulaire, connue de I'homme du 
metier, pour les filtres orthogonaux. 

Pour realiser des filtres en ondelettes avec un gain DC unite et un 
gain de Nyquist de 2, c'est-a-dire avec une normalisation dite (1,2), aucune 
20 etape specifique de normalisation du type Xj = G.x i( ou G est une valeur reelle, 
n'est necessaire. En revanche, si une normalisation differente est requise, telle 
qu'une normalisation (Gi, 2G 2 ) avec un gain DC de Gi et un gain de Nyquist de 
2G2, une etape de normalisation inverse est necessaire : x 2n =( 1/ G 1 ).x 2n et 

x 2n+1 =(1/G 2 )-X 2 n+1- 
25 Comme la montre la figure 6, on procede ensuite a une iteration 

d'etapes de lifting passe-bas WSA (etape 602) et passe-haut WSA (etape 608), 
soit k max -1 etapes de lifting WSA au total, la valeur du compteur k decroissant 

de k max -1 a 1. 
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Pour cela, on decremente d'une unite la valeur du compteur k apres 
chaque passage par I'etape 602 (voir etape 604 sur la figure 6) et apres chaque 
passage par I'etape 608 (voir etape 610 sur la figure 6). 

On reitere les etapes 602, 604, 608 et 610 jusqu'a ce que la valeur 
5 de k soit egale a 1 (voir tests 606 et 612 sur la valeur de k, respectivement 
effectues apres les etapes de decrementation 604 et 610). 

Les etapes de lifting passe-bas WSA 602 sont definies par : 

p-2k-1-1 

Vn,X 2 n = X 2n "R I>2k-1,j-(*2n+1-2j " *2n+1+2j) ( 27 ) 

I i=o J 
et les etapes de lifting passe-haut WSA 608 sont definies par : 

fL 2k -1 >| 
10 Vn,x 2n+1 = x 2n+1 -R l>2k,j-( x 2n-2j -x 2n+2j ) (28) 

U=° J 
Ces etapes sont connues en soi. 

Lorsque le resultat du test 606 ou 612 est positif, c'est-a-dire que k = 
1 , I'alternance d'etapes de lifting passe-bas WSA et passe-haut WSA prend fin 
et est suivie d'une etape de lifting passe-bas RSC 614, definie par : 
15 Vn,x 2n =x 2n -R(y.x 2n+1 ) (29) 

Dans un mode de realisation prefere, le coefficient y est calcule en 
fonction du coefficient de ponderation utilise a I'etape de lifting passe-haut 
suivante : 

Y = -1/(2Po,L 0 -l) < 30 > 
20 Lors de I'etape 614, on initialise aussi le parametre m } a la valeur 1 et 

un compteur j a la valeur L 0 -1. 

A Tissue de I'etape de lifting passe-bas 614, on effectue un test 616 
pour determiner si la valeur du compteur j est egale ou non a 0. 

Tant que le resultat du test 616 est negatif, on itere une etape 618 
25 qui consiste essentiellement en une sequence d'operations de modification et 
de ponderation des echantillons de rang pair ou impair du signal x n . 

En effet, lors de I'etape 618, on effectue une succession de L 0 
couples d'etapes de lifting, correspondant a des valeurs du compteur j 
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comprises entre L 0 -1 et 0, ou la premiere etape de lifting est une etape de lifting 
passe-haut HSA : 

Vn,x 2n+1 = x 2n+ i -R(Po,j*(X2n -*2n + 2)) (31) 

et ou la seconde etape de lifting est une etape de lifting passe-bas a coefficient 
5 unique : 

Vn,x 2n = x 2n -R(aoj-X2n+mj ) ( 32 ) 

Lorsque le resultat du test 616 est positif, c'est-a-dire que j = 0, on 
effectue la derniere etape de lifting 620, qui est I'inverse de la premiere etape 
de lifting passe-haut effectuee dans la transformation directe. Dans un mode 
10 prefere de realisation, cette etape de lifting est du type, soit LSC, soit RSC, en 

fonction de la valeur du parametre m 0 = (-1) L ° , et est definie par : 

Vn,x 2n+1 = x 2n+1 +R(x 2n+1 _ mQ ) (33) 

La figure 7 illustre une procedure permettant de determiner les 
parametres de lifting pour tout couple de filtres en ondelettes HSS/HSA donne, 

15 defini par les equations (4), (6) et (7) figurant en introduction. Cette procedure 
s'applique aussi aux filtres orthogonaux ainsi qu'aux filtres quelconques, sans 
symetrie particuliere, de support egal. Les procedures "UNLIFT" apparaissant 
ci-dessous sont detaillees plus loin. 

Comme le montre la figure 7, une premiere etape 700 consiste a 

20 initialiser une variable entiere kaOet une variable n 0 a -1 et a memoriser dans 
une variable Ij, i <e {0,1}, la longueur, c'est-a-dire le nombre de coefficients, des 
filtres passe-bas et passe-haut intermediates H 0 et Hi, respectivement (c'est-a- 
dire le nombre de coefficient des filtres dits equivalents). 

Puis une etape 702 consiste a tester si lo < h- 

25 Si le resultat du test 702 est positif, on effectue une etape 704 de 

determination des coefficients de f etape de lifting passe-haut WSA 518 de la 
figure 5 en appliquant une procedure UNLIFT_H1 (voir figure 8B), a Tissue de 
laquelle on obtient le nombre de coefficients de lifting L k , les L k coefficients a k j 
pour 0 < i < L k et le nouveau filtre passe-haut intermediate Hi. 
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A Tetape 704, on incremente aussi d'une unite la valeur de k, on met 
a jour la longueur lj du filtre intermediaire et on memorise la valeur 1 dans la 
variable no. 

Si le resultat du test 702 est negatif, on effectue une etape 706 
5 consistant a tester si l 0 > h. 

Si le resultat du test 706 est positif, on effectue .une etape 708 de 
determination des coefficients de I'etape de lifting passe-bas WSA 512 de la 
figure 5 en appliquant une procedure UNLIFTH0 (voir figure 8A), a Tissue de 
laquelle on obtient le nombre de coefficients de lifting L k , les L k coefficients a k j 
10 pour 0 < i < L k et le nouveau filtre passe-bas intermediaire Ho. 

A Tetape 708, on incremente aussi d'une unite la valeur de k, on met 
a jour la longueur lj du filtre intermediaire et on memorise la valeur 0 dans la j 
variable n 0 . 

Si le resultat du test 706 est negatif, on effectue une etape 710 
1 5 consistant a tester si la variable n 0 est nulle. 

Si le resultat du test 710 est positif, on effectue une etape 712 
consistant a memoriser la valeur 0 dans la variable L k et a incrementer d'une 
unite la valeur de k. 

A Tissue de Tetape 712, ou bien si le resultat du test 710 et negatif, 
20 on effectue une etape 714 de determination des coefficients de Tetape de lifting 
passe-bas a coefficient unique incluse dans Tetape 506 de la figure 5 en 
appliquant une procedure UNLIFT_EQ (voir figure 8C), a Tissue de laquelle on 
obtient le nombre de coefficients de lifting L 0 et les L 0 coefficients a 0 ,i pour 
0 < i < L 0 . 

25 A Tetape 714, on incremente aussi d'une unite la valeur de k. 

A I'issue de Tetape 714, on effectue une etape 716 de 
rearrangement, consistant a mettre a jour les coefficients de lifting obtenus en 
les reordonnant, comme suit : 

(Lo.Li, ...,Lk-i) = (Lk-i,Lk-2.-..,Lo) 
30 (a 0 ,i,ai i j,...,a k -i i i) = (a k -ij,a k - 2( i,... ) ao l j) 

Lors de Tetape 716, on memorise aussi la valeur k max dans la 

variable k. 
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En resume : 

- lorsque l 0 > h, les coefficients de lifting sont obtenus par la 
procedure UNLIFT_HO(H 0 ,H 1 ), 

- lorsque l 0 < h, les coefficients de lifting sont obtenus par la 
5 procedure UNLIFT^HIfHo.HO et 

- lorsque l 0 = li, les coefficients de lifting sont obtenus par la 
procedure UNLIFT_EQ(H 0l Hi). 

On decrit maintenant ces trois procedures plus en detail. 

Comme le montre la figure 8A, qui illustre la procedure UNLIFT_H0, 
10 une etape d'initialisation 800A consiste a memoriser dans la variable l 0 la moitie 
de la longueur de H 0 et a memoriser dans la variable li la moitie de la longueur 
de Hv Lors de I'etape 800A, on initialise aussi une variable i et la variable L k a 
la valeur 0. 

Puis une etape 802A consiste a tester si l 0 > h. 
15 Si le resultat du test 802A est negatif, la procedure UNLIFTJHO est 

terminee. 

Si le resultat du test 802A est positif, une etape 804A consiste a 
calculer a kt i et H 0 (z), comme suit : 

a k,i = " H 0,l 0 /H Ui 

H 0 (z) = H 0 (z) - a kJ .(z"'o + ' 1 - z'o "'i (z) 

20 En fin d'etape 804A, on decremente de deux unites la valeur de l 0 et 

on incremente d'une unite la valeur de i et la valeur de L k . 

Puis on reitere cette etape 804A tant que l 0 > li. 
Comme le montre la figure 8B, qui illustre la procedure UNLIFTJH1, 
une etape d'initialisation 800B identique a I'etape 800A de la figure 8A est 
25 d'abord effectuee. 

Puis une etape 802B consiste a tester si U > l 0 . 

Si le resultat du test 802B est negatif, la procedure UNLIFT H1 est 

terminee. 

Si le resultat du test 802B est positif, une etape 804B consiste a 
30 calculer otk.j et Hi(z), comme suit : 
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<*k,i = " H 1.h /H 0J 0 

H 1 (z)= s H 1 (z)-a kii .(z- | i +, o -z'i-'oj.HqCz) 

En fin d'etape 804B, on decremente de deux unites la valeur de h et 
on incremente d'une unite la valeur de i et la valeur de L k . 

Puis on reitere cette etape 804B tant que U > lo- 
5 Comme le montre la figure 8C, qui illustre la procedure UNLIFT^EQ, 

une etape d'initialisation 800C consiste a memoriser dans la variable lo la moitie 
de la longueur de Ho, a memoriser dans la variable U la valeur de lo- Lors de 
I'etape 800C, on initialise aussi la variable L 0 a la valeur l 0 -1 et une variable i a 
la valeur L 0 -1. 

10 L'etape suivante 802C consiste a calculer H 0 (z) et le parametre y, 

comme suit : 

y = H oj 0 /H L'i 

H 0 (z) = z.[H 0 (z)- Y .H 1 (z)] 

En fin d'etape 802C, on decremente d'une unite la valeur de l 0 et on 
initialise une variable q a la valeur 0. 
15 Puis une etape 804C consiste a tester si l 0 > 0. 

Si le resultat du test 804C est negatif, la procedure UNLIFT_EQ est 

terminee. 

Si le resultat du test 804C est positif, une etape 806C consiste a 
tester si q > 0. 

20 Si le resultat du test 806C est positif, I'etape suivante 808C consiste 

a calculer p 0 ,i, H^z), a 0 ,i et H 0 (z), comme suit : 

Po.i =H 1,-I 1 +1 /H 0-I 0 +1 

H 1 (z) = z- 1 .[H 1 (z)-p 0ii .(z- 1 -z).H 0 (z)] 

a 0,i ^- H 0-l 0 +1 /H 1-l-,+1 

H 0 (z) = z- 1 .[Ho(z)-a 0i i.H 1 (z)] 

Au cours de I'etape 808C, apres calcul de Hi(z), on decremente 
d'une unite la valeur de U ©t, apres calcul de H 0 (z), on decremente d'une unite 
25 la valeur de l 0 et la valeur de i et on memorise la valeur 1 dans la variable q. 
Puis on reitere cette etape 808C tant que l 0 > 0. 
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Si le resultat du test 806C est negatif, I'etape suivante 81 OC consiste 
a calculer p 0 ,j, H^z), a 0 ,i et H 0 (z), comme suit : 

Po.i =- h 1,I 1 /H 0,l 0 

H 1 (z) = z.[H 1 ( 2 )-(3 0ii .(z- 1 -z).H 0 (z)] 

a O,i = H 0,I 0 ""h^ 

H o (z) = z.[H 0 (z)-a Ot] Hi(z)) 

Au cours de I'etape 81 OC, apres calcul de H^z), on decremente 
d'une unite la valeur de U et, apres calcul de H 0 (z), on decremente d'une unite 
la valeur de l 0 et la valeur de i et on memorise la valeur 0 dans la variable q. 

Puis on reitere cette etape 810C tant que l 0 > 0. 

Les figures 9 et 10 illustrent sous forme schematique un bloc de 
transformation elementaire compris dans un dispositif de filtrage conforme a la 
presente invention, dans deux modes particuliers de realisation. La figure 9 
illustre plus particulierement un bloc de transformation par lifting directe et la 
figure 10 illustre plus particulierement un bloc de transformation par lifting 
inverse. 

Le dispositif de filtrage 900 schematise sur la figure 9 est equivalent 
a celui du bloc BE illustre sur la figure 4. 

II comporte une entree E 2 o a laquelle est applique le signal a 
transformer. Le signal a transformer comporte une suite d'echantillons {Xj}, ou i 
est un indice de rang d'echantillon. 

L'entree E 2 o est reliee a un premier decimateur par deux D 2 o qui 
delivre les echantillons {x 2 i} de rang pair. 

L'entree E 20 est egalement reliee a une avance AV 21 suivie d'un 
second decimateur D 21j qui delivre les echantillons {x 2i+1 } de rang impair. 

Le premier decimateur D 20 est relie a un premier filtre A 0 qui filtre les 
echantillons de rang pair en leur appliquant un coefficient de ponderation p 0 ,j et 
une fonction d'approximation comme decrit dans I'equation (18) et les fournit a 
un additionneur AD 2 o- 

L'ensemble forme par le filtre A 0 et I'additionneur AD 2 o est un module 
de lifting passe-haut, note FHP sur le dessin. 



34 

La sortie du second decimateur D 2 i est egalement reliee a 
I'additionneur AD 2 o. 

En consequence, I'additionneur AD20 delivre des echantillons 
intermediates y 2i+ i calcules selon I'equation (18). 
5 La sortie de I'additionneur AD 2 o est egalement reliee a un second 

filtre A1 qui filtre les echantillons de rang impair en leur appliquant un coefficient 
de ponderation a 0 j et une fonction d'approximation comme decrit dans 
I'equation (17) et les fournit a un additionneur AD 2 l 

L'ensemble forme par le filtre et I'additionneur AD21 est un module 
1 0 de lifting passe-bas, note FLP sur le dessin. 

La sortie du premier decimateur D 20 est egalement reliee a 
I'additionneur AD21. 

La sortie de I'additionneur AD21 est la premiere sortie S 20 du bloc de 
transformation, qui delivre un signal numerique comportant des echantillons 
1 5 basse frequence y 2 \. 

L'additionneur AD 2 o est relie a une seconde sortie S 2 i du bloc 
elementaire, qui delivre un signal numerique comportant des echantillons haute 
frequence y 2 i+i- 

Le dispositif de filtrage 1000 schematise sur la figure 10 effectue la 
20 transformation inverse de cede effectuee par le dispositif de filtrage 900 
schematise sur la figure 9. 

Ce bloc de transformation comporte une premiere entree E 30 a 
laquelle est applique un premier signal a transformer et une seconde entree E31 
a laquelle est applique un second signal a transformer 
25 Les signaux a transformer comportent ici les echantillons obtenus 

apres filtrage d'analyse d'un signal numerique par le bloc d'analyse de la figure 
9. Ces echantillons ont eventuellement ete modifies par un autre traitement 
entre I'analyse et la synthese. 

Plus precisement, le premier signal a transformer comporte les 
30 echantillons basse frequence {y 2 i} et le second signal a transformer comporte 
les echantillons haute frequence{y 2j+ i}. 
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Le bloc de synthese a une structure analogue a celle du bloc 
d'analyse, et peut en etre deduit simplement. En particulier, le bloc de synthese 
utilise les memes filtres A 0 et A-i que le bloc d'analyse. 

L'entree E 30 est reliee a un soustracteur SO 30 . 
5 L'entree E 31 est reliee au filtre Ai, lui-meme relie au soustracteur 

SO 30 . 

L'ensemble forme par le filtre et le soustracteur SO 30 est un 
module de lifting passe-bas, note ILP sur le dessin. 

La sortie du soustracteur SO 30 delivre un signal numerique 
10 comportant des echantillons intermediates y 2 \ calcules selon Pequation (32). 

La sortie du soustracteur SO 30 est reliee au filtre A 0j lui-meme relie a 
un soustracteur S0 3 i auquel est egalement reliee l'entree E 31 . 

L'ensemble forme par le filtre A 0 et le soustracteur S0 3 i est un 
module de lifting passe-haut, note IHP sur le dessin. 
15 La sortie du soustracteur S0 3 i delivre un signal numerique 

comportant des echantillons reconstruits de rang impair x 2 i+i. 

La sortie du soustracteur SO 30 delivre un signal numerique 
comportant des echantillons reconstruits de rang pair x 2 \. 

La sortie du soustracteur SO 30 est reliee a un interpolateur par deux 
20 IN 30 et la sortie du soustracteur S0 3 i est reliee a un interpolateur par deux IN 31 
lui-meme relie a un retard R 31 . L'interpolateur IN 30 et le retard R 3 i sont relies a 
un additionneur AD 32 qui delivre en sortie S 30 le signal comportant les 
echantillons Xj reconstruits. 

On donne maintenant un exemple, nullement limitatif, de 
25 transformation directe (en anglais "forward') et inverse par lifting pour la 
realisation de filtres en ondelettes conformement a I'invention. 

La paire de filtres utilisee dans cet exemple est la suivante : 

H 0 (z)= a - Z ~ 1+1 + Z + a - z2 (35) 
U 2.(a + 1) V 7 

HiM- "" 1 ;, 1 -^'" 2 (36) 
2.(a-1) 

30 Les parametres de lifting correspondant a cette paire de filtres sont : 
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k = 1 , L 0 = 1 et cto.o = a- (37) 
On choisit egalement une fonction d'approximation R identique pour 

toutes les etapes de lifting et definie par R(x) = E(x+1/2) ou E designe la partie 

entiere. 

5 La longueur du signal est choisie egale a 6 : i 0 = 0 et h = 6. 

Pour ce qui concerne les equations de lifting directes, toutes les 
valeurs du signal de sortie Y sont initialisees aux valeurs du signal d'entree X. 
Le signal Y est ensuite etendu d'un echantillon a chaque limite : y_i = y 0 et y6 = 

10 Les quatre operations de lifting suivantes sont effectuees. Tout 

d'abord, les valeurs des echantillons (y.i, y 1f y 3) y 5 ) sont calculees a partir de : 

y 2n+ i = y 2 n + i-E(y 2 n + 2+1/2) pour n = -1, 0, 1, 2 (38). 
Ensuite, les valeurs des echantillons (y 0 , y 2l y4, y6> sont calculees a 

partir de : 

15 y 2n = y 2 n+E(a 0 , 0 -y 2 n-i+1/2) pour n = 0, 1, 2, 3 (39) 

ou 

ao.o = cx (40) 
Puis les valeurs des echantillons (yi, y 3 , y 5 ) sont calculees a partir 

de : 

20 y 2n+ i = y 2 n + i+E[p 0l0 .(y 2 n-y 2 n +2 )+1/2)] pour n = 0, 1, 2 (41) 

ou 

Po,o = -1/(1-2a) (42) 
Enfin, les valeurs des echantillons (y 0l y 2 , y4) sont calculees a partir 

de : 

25 y 2n = y 2 n+E(y.y 2n+ i+1/2) pour n = 0, 1 , 2 (43) 

ou 

Y = -1/(2p 0 ,o) (44) 
Pour ce qui concerne les equations de lifting inverses, toutes les 
valeurs du signal de sortie X sont initialisees aux valeurs du signal d'entree Y. 
30 Le signal X est ensuite etendu de deux echantillons a chaque limite : x. 2 = xo, x_i 
= -xi et x 6 = Xa et x 7 = -x 5 . 
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Les quatre operations de lifting inverse suivantes sont effectuees. 
Tout d'aborcl, les valeurs des echantillons (x_ 2l x 0 , x 2> x 4j x 6 ) sont calculees a 
partir de : 

x 2n = x 2 n-E(y.x 2 n + i+1/2) pour n = -1 , 0, 1 , 2, 3 (45) 
5 Ensuite, les valeurs des echantillons (x. 1( x 1t x 3) x 5 ) sont calculees a 

partir de : 

x 2n +i = x 2 n + i-E[po ( o.(x 2 n-x 2 n + 2) + 1/2)] pour n = -1 , 0, 1 , 2 (46) 
Puis les valeurs des echantillons (x 0 , x 2 , x 4 , x 6 ) sont calculees a partir 

de : 

10 x 2n = x 2n -E(a 0 ,o.x 2n -i + 1/2) pour n = 0, 1, 2, 3 (47) 

Enfin, les valeurs des echantillons (x 1( x 3 , x 5 ) sont calculees a partir 

de : 

x 2n +i = x 2n+ i+E(x 2n+2 +1/2) pour n = 0, 1 , 2 (48) 
Le champ d'application de la presente invention ne se limite pas aux 
15 filtres en ondelettes HSS/HSA de support egal, mais s'etend beaucoup plus 
largement, d'une part, aux filtres en ondelettes orthogonaux et, d'autre part, a 
des filtres en ondelettes quelconques de support egal. 

En ce qui concerne I'application de ('invention a la realisation par 
lifting de filtres en ondelettes orthogonaux, la seule difference par rapport a ce 
20 qui a ete decrit precedemment reside dans la valeur des parametres de lifting 
ccoj, po.i, y et dans les relations entre ces parametres, afin de garantir qu'ils 
realisent indifferemment des filtres HSS/HSA ou des filtres orthogonaux. II est a 
noter que la methode de calcul des parametres de lifting, decrite plus haut en 
liaison avec les figures 7 et 8, peut etre utilisee pour engendrer tous les 
25 parametres de lifting a 0( i, poj, y pour les filtres HSS/HSA aussi bien que pour les 
filtres orthogonaux. 

Tous les filtres orthogonaux sont de support egal et de longueur 
paire mais ne sont pas symetriques, sauf dans un cas trivial : le filtre de Haar, 
defini par H 0 (z) = (1+z)/2 et Hi(z) = -1+z. 
30 Les filtres orthogonaux de longueur 4 sont caracterises par un 

parametre 8 representant un degre de liberie : 
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-1 2 
cos 0 (sin G + cos 0).z + sin 0.(sin G + cos G) + sin 0(sin 0 - cos G).z + cos O.(cos G - sin G).z 

H 0 ( Z ) = (49) 

2 

H-| (z) = cosG.(cosG - sinO).z" 1 - sinG.(sinG - cosG) + sinG.(sinO + cosG).z - cosG.(sinG + cosG).z 2 (50) 

De fa<?on plus generale, les filtres orthogonaux satisfont la relation : 

H 1 (z) = 2z.H 0 (-z- 1 ) (51) 

5 Dans I'exemple non limitatif donne ici, on choisit, de meme que 

precedemment, une fonction d'approximation R identique pour toutes les 

etapes de lifting et definie par R(x) = x+1/2. 

La longueur du signal est choisie egale a 6 : i 0 = 0 et h = 6. 

La realisation par lifting d'un tel filtre se fait en utilisant les memes 

10 etapes que dans I'exemple precedent (excepte pour ce qui concerne la 

procedure d'extension, decrite plus loin) a ceci pres que les parametres sont ici 

choisis de la fagon suivante : 

a . a cos 9 _ sin 6 -cos 9 /cox 

k= 1, L 0 = 1, a 00 =— — Po.o = cos(29), y = —— — (52) 

■ sin9 + cos9 2(sin9 + cos0) 

La procedure de factorisation lifting donnee plus haut est similaire 
15 pour les filtres orthogonaux. La seule difference reside dans la relation entre les 
coefficients aoj, poj. Y- 

Comme on I'a vu pour les filtres HSS/HSA, les parametres otoj sont 
fonction des coefficients des filtres de la transformation en ondelettes a realiser. 
II en est de meme pour les filtres orthogonaux. La methode de calcul des 
20 parametres de lifting decrite plus haut est inchangee. 

Comme on I'a vu pour les filtres HSS/HSA, les parametres p 0l i sont 
fonction des parametres a 0( i et sont tels qu'ils garantissent que le filtre passe- 
haut Hi(z) equivalent de chaque etape de lifting passe-haut soit HSA. La 
situation est analogue pour les filtres orthogonaux : les parametres p 0 j sont 
25 fonction des parametres a 0| i et sont tels qu'ils garantissent que le filtre passe- 
haut H^z) equivalent de chaque etape de lifting passe-haut soit tel que 

H 0 (z) = -z.H 1 (-z" 1 )/2 et H^z) definissent une transformation en ondelettes 
orthogonale. La methode de calcul des parametres de lifting decrite plus haut 
peut aussi etre utilisee pour calculer les parametres poj. 
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Enfin, comme on I'a vu pour les filtres HSS/HSA, le parametre y est 
fonction des parametres Po,l 0 -i et est tel ^ u il garantisse que le filtre passe- 
bas H 0 (z) equivalent de cette etape de lifting passe-bas soit HSS. La situation 
est analogue pour les filtres orthogonaux : le parametre y est fonction des 
5 parametres a 0 ,i et p 0 ,i et est tel qu'il garantisse que le filtre passe-bas H 0 (z) 
equivalent de cette etape de lifting passe-bas soit tel que 

H 0 (z) = -z.H 1 (-z" 1 )/2. La methode de calcul des parametres de lifting decrite 

plus haut peut aussi etre utilisee pour calculer le parametre y. 

Pour la transformation par lifting directe, toutes les valeurs du signal 
10 de sortie Y sont initialisees aux valeurs du signal d'entree X. Le signal Y est 
ensuite etendu d'un echantillon a chaque limite : dans I'exemple fourni ici, on 
utilise une procedure d'extension connue de I'homme du metier sous le nom 
d'extension circulaire : y_i = y 5 et y 6 = yo- 

Toutes les etapes ulterieures sont identiques a celles de I'exemple 
15 precedent, en utilisant les valeurs donnees a I'equation (52). 

Pour la transformation par lifting inverse, toutes les valeurs du signal 
de sortie X sont initialisees aux valeurs du signal d'entree Y. Le signal X est 
ensuite etendu de deux echantillons a chaque limite : x. 2 = x 4l x_i = x 5l x 6 = x 0 et 

x 7 = Xi. 

20 Toutes les etapes ulterieures sont identiques a celles de I'exemple 

precedent, en utilisant les valeurs donnees a ('equation (52). 

En ce qui concerne I'application de 1'invention a la realisation par 
lifting de filtres en ondelettes quelconques de support egal, les seules 
differences par rapport a ce qui a ete decrit precedemment pour les filtres 

25 HSS/HSA et les filtres orthogonaux sont les suivantes : 

- il n'y a pas de relation de dependance entre les parametres a 0 ,i, 
po.i, y ; il n'est done pas necessaire de calculer ces valeurs dans la RAM ; 

- la premiere etape de lifting passe-haut dans la realisation par lifting 
directe peut etre quelconque, e'est-a-dire : 



30 y 2n +1 =Y2n + 1- R 



X § i-y2n + 2i 



V i 



(53) 
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ou 5j est un coefficient de ponderation. 

II est a noter que par exemple, si on utilise dans la transformation 
directe, en lieu et place de I'equation (16), I'equation suivante : 




5 alors, dans la transformation inverse, la derniere etape de lifting 620 est 
exactement son inverse, soit, dans le cas ci-dessus : 



realisation par lifting d'un couple de filtres en ondelettes quelconques de 
10 support egal peuvent etre obtenus par I'intermediaire de la methode de calcul 
des parametres de lifting decrite plus haut en liaison avec les figures 7 et 8. Les 
coefficients de lifting 6j sont alors les coefficients d'indice pair du filtre H^z) 
obtenu a partir de la formule (L k , a k(it Hi, H 0 ) = UNLIFT_EQ(H 0) H 1 ) de I'etape 
714 illustree sur la figure 7. 




(34) 



Tous les coefficients de lifting a 0 ,i, pi 



oj, y necessaires pour la 




REVENDICATIONS 



1. Procede de filtrage adapte a transformer un signal numerique 
d'entree (x n ) en un ou plusieurs signaux numeriques de sortie (y n ) comportant 

5 des echantillons de rang pair (y 2n ) et des echantillons de rang impair (y 2 n+i), 
ledit procede comportant au moins une iteration (506) qui contient 

- une operation de modification d'echantillons de rang pair (y 2n ) par 
une fonction (R) d'echantillons de rang impair ponderes (a 0 j.y2n+mj ). 

- une operation de modification d'echantillons de rang impair (y 2n +i) 
10 par une fonction (R) d'echantillons de rang pair ponderes (Po,j (y2n -Y2n+2)). 

lesdits echantillons ponderes etant obtenus par au moins une 
operation de ponderation, 

ledit procede etant caracterise en ce qu'au moins une desdites operations de 
ponderation est appliquee a la difference entre deux echantillons de rang pair 
15 consecutifs. 

2. Procede de filtrage selon la revendication 1, caracterise en ce que 
ladite operation de modification d'echantillons de rang impair (y 2n +i) est 
effectuee a la suite de ladite operation de modification d'echantillons de rang 
pair (y 2n ). 

20 3. Procede de filtrage selon la revendication 1 ou 2, caracterise en ce 

que ladite iteration (506) consiste notamment a : 

- ponderer, par un premier coefficient de ponderation (a 0 j), au 

moins un echantillon de rang impair (y2n+mj ) voisin d'un echantillon de rang 

pair en cours de modification, de fagon a obtenir un echantillon de rang impair 
25 pondere (a 0i j.y2n+mj ) , 

- modifier au moins un echantillon de rang pair (y 2n ) a partir d'au 
moins un echantillon de rang impair pondere (a 0i j.y2n+mj 

- ponderer, par un deuxieme coefficient de ponderation (Pqj), des 
echantillons de rang pair (y 2n -Y2n-f2) voisins d'un echantillon de rang impair 




en cours de modification, de facpon a obtenir des echantillons de rang pair 
ponderes (P 0i j(y2n -y2n + 2» - et 

- modifier au moins un echantillon de rang impair (y2n+i) a partir d'au 
moins un echantillon de rang pair pondere (Po t j.(y2n - Y2n+2 )) • 

5 4. Procede de filtrage selon la revendication precedente, caracterise 

en ce que le deuxieme coefficient de ponderation *(p 0 j) est fonction du premier 

coefficient de ponderation (a 0 j). 

5. Procede de filtrage selon la revendication precedente, caracterise 
en ce que le deuxieme coefficient de ponderation (Poj) depend du premier 

10 coefficient de ponderation (cxqj) comme suit : 



j ^ 



1-22> 0il 
i=0 



ou ao,i designe le premier coefficient de ponderation, p 0 j designe le deuxieme 
coefficient de ponderation, i et j sont des entiers et rrij est une valeur definie par 

la recurrence mo = (-1) L ° et mj = -nrij.i, L 0 etant un entier predetermine. 

15 6. Procede de filtrage selon Tune quelconque des revendications 

precedentes, caracterise en ce que, a chaque iteration, I'echantillon de rang 
impair (y2n+mj ) voisin de I'echantillon pair en cours de modification est 

alternativement I'echantillon de rang immediatement inferieur (y2n-i) ou 

immediatement superieur (y2n+i). 
20 7. Procede de filtrage selon Tune quelconque des revendications 

precedentes, caracterise en ce qu'il comporte, a Tissue de ladite iteration (506), 

une etape de filtrage supplemental (508) incluant une operation de 

ponderation par un troisieme coefficient de ponderation (y). 

8. Procede de filtrage selon la revendication precedente, caracterise 
25 en ce que le troisieme coefficient de ponderation (y) est fonction du coefficient 

de ponderation utilise a I'etape precedente, comme suit : 

Y = -1/(2Po,L 0 -l) 
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ou y designe le troisieme coefficient de ponderation, L 0 est un parametre 
predetermine et Po,i_o-i designe le coefficient de ponderation utilise a I'etape 
precedente. 

9. Procede de filtrage selon Tune quelconque des revendications 
5 precedentes, caracterise en ce que le signal numerique d'entree (x n ) represente 

une image. 

10. Procede de filtrage adapte a transformer un ou plusieurs signaux 
numeriques d'entree (y n ) en un signal numerique de sortie (x n ), lesdits signaux 
d'entree comportant des echantillons de rang pair (y 2n ) et des echantillons de 

10 rang impair (y 2 n+i), ledit procede comportant au moins une iteration (618) qui 
contient 

- une operation de modification d'echantillons de rang impair (x 2n +i) 
par une fonction (R) d'echantillons de rang pair ponderes (Poj.(x2n " x 2n+2)) . 

- une operation de modification d'echantillons de rang pair (x 2n ) par 
15 une fonction (R) d'echantillons de rang impair ponderes (a 0 i j.(x 2 n +m j ))» 

lesdits echantillons ponderes etant obtenus par au moins une 
operation de ponderation, 

ledit procede etant caracterise en ce qu'au moins une desdites operations de 
ponderation est appliquee a la difference entre deux echantillons de rang pair 
20 consecutifs. 

11. Procede de filtrage selon la revendication precedente, 
caracterise en ce que ladite operation de modification d'echantillons de rang 
pair (x 2n ) est effectuee a la suite de ladite operation de modification 
d'echantillons de rang impair (x 2n +i). 

25 12. Procede de filtrage selon la revendication 10 ou 11, caracterise 

en ce que ladite iteration (618) consiste notamment a : 

- ponderer, par un quatrieme coefficient de ponderation (Pop, des 

echantillons de rang pair (x 2n ~ x 2n+2) voisins d'un echantillon impair en cours 
de modification, de fa?on a obtenir des echantillons de rang pair ponderes 
30 (Po.r(><2n-x 2n+2 )) ( 




- modifier au moins un echantillon de rang impair (x 2n +i) a partir d'au 
moins un echantillon de rang pair pondere (Po, j.(*2n - x 2n+2 )) . 

- ponderer, par un cinquieme coefficient de ponderation (txQj), au 

moins un echantillon de rang impair (x 2 n+mj ) voisin d'un echantillon pair en 

5 cours de modification, de fa?on a obtenir un echantillon de rang impair pondere 
(aoj-^n+mj ) et 

- modifier au moins un echantillon de rang pair (x 2n ) a partir d'au 
moins un echantillon de rang impair pondere (a 0i j.x 2n+mj ). 

13. Procede de filtrage selon la revendication precedente, 
10 caracterise en ce que le quatrieme coefficient de ponderation (Poj) est fonction 

du cinquieme coefficient de ponderation (otQj). 

14. Procede de filtrage selon la revendication precedente, 
caracterise en ce que le quatrieme coefficient de ponderation (Poj) depend du 

cinquieme coefficient de ponderation (a 0( j) comme suit : 



f 



15 p oj = mj / 



j 

1-2Za 0ii 

i=0 j 



ou aoj designe le cinquieme coefficient de ponderation, Po.j designe 
le quatrieme coefficient de ponderation, i et j sont des entiers et rrij est une 

valeur definie par la recurrence m 0 =(-1) L ° et mj = -mj.i, L 0 etant un entier 

predetermine. 

20 15. Procede de filtrage selon Tune quelconque des revendications 10 

a 14, caracterise en ce que, a chaque iteration, I'echantillon de rang impair 
( x 2n+mj ) voisin de I'echantillon pair en cours de modification est 

alternativement Techantillon de rang immediatement inferieur (x 2n -i) ou 
immediatement superieur (x 2n +i). 
25 16. Procede de filtrage selon Tune quelconque des revendications 10 

a 15, caracterise en ce qu'il comporte, prealablement a ladite iteration (618), 



45 



une etape de filtrage supplemental (614) incluant une operation de 
ponderation par un sixieme coefficient de ponderation (y). 

17. Procede de filtrage selon la revendication precedente, 
caracterise en ce que le sixieme coefficient de ponderation (y) est fonction du 
5 coefficient de ponderation utilise a I'etape suivante, comme suit : 

Y = -1/(2P 0t Lo-l) 

ou y designe le sixieme coefficient de ponderation, L 0 est un parametre 
predetermine et Po,l 0 -i d6si 9 ne le coefficient de ponderation utilise a I'etape 
suivante. 

10 18. Procede de filtrage selon Tune quelconque des revendications 10 

a 17, caracterise en ce que le signal numerique de sortie (x n ) represente une 
image. 

19. Procede de filtrage selon Tune quelconque des revendications 
precedentes, caracterise en ce que lesdites operations de modification 

15 consistent a appliquer une fonction d'approximation (R). 

20. Procede de filtrage selon la revendication precedente, 
caracterise en ce que la fonction d'approximation (R) est la fonction identite. 

21. Procede de filtrage selon la revendication 19, caracterise en ce 
que la fonction d'approximation (R) est une fonction d'une variable reelle qui 

20 fournit rentier le plus proche de la variable. 

22. Procede de filtrage selon la revendication 19, caracterise en ce 
que la fonction d'approximation (R) est une fonction d'une variable reelle qui 
fournit le premier entier inferieur a la variable. 

23. Procede de filtrage selon la revendication 19, caracterise en ce 
25 que la fonction d'approximation (R) est une fonction d'une variable reelle qui 

fournit le premier entier superieur a la variable. 

24. Procede de filtrage selon la revendication 19, caracterise en ce 
que la fonction d'approximation (R) est une fonction d'une variable decomposee 
en sous-variables dont la somme est egale a la variable, qui fournit une somme 

30 de valeurs approchees des sous-variables, chacune des valeurs approchees 
des sous-variables etant, soit une fonction d'une variable reelle qui fournit 




rentier le plus proche de la variable, soit une fonction d'une variable reelle qui 
fournit le premier entier inferieur a la variable, soit une fonction d'une variable 
reelle qui fournit le premier entier superieur a la variable. 

25. Dispositif de traitement de signal (10), caracterise en ce qu'il 
5 comporte des moyens adaptes a mettre en oeuvre un procede de filtrage selon 

Tune quelconque des revendications precedentes. 

26. Dispositif de filtrage numerique adapte a transformer un signal 
numerique d'entree (x n ) en un ou plusieurs signaux numeriques de sortie (y n ) 
comportant des echantillons de rang pair (y 2n ) et des echantillons de rang 

10 impair (y 2 n+i), ledit dispositif de filtrage comportant 

- au moins un moyen de ponderation, 

- des moyens pour modifier des echantillons de rang pair (y 2n ) par 
une fonction (R) d'echantillons de rang impair ponderes (aoj-V^n+mj ) . 

- des moyens pour modifier des echantillons de rang impair (y2n+i) 
15 par une fonction d'echantillons de rang pair ponderes (Poj-(y2n ~ y2n+2 )) - 

lesdits echantillons ponderes etant fournis par ledit moyen de 
ponderation, lesdits moyens de modification fonctionnant de fagon iterative, de 
fa?on a modifier au moins une fois des echantillons de rang pair (y 2n ) puis au 
moins une fois des echantillons de rang impair (y 2 n+i), 
20 ledit dispositif de filtrage etant caracterise en ce qu'au moins un desdits moyens 
de ponderation recpoit en entree la difference entre deux echantillons de rang 
pair consecutifs. 

27. Dispositif de filtrage selon la revendication precedente, 
caracterise en ce que lesdits moyens pour modifier des echantillons de rang 

25 impair (y 2n +i) sont disposes en aval desdits des moyens pour modifier des 
echantillons de rang pair (y 2n ). 

28. Dispositif de filtrage selon la revendication 26 ou 27, caracterise 
en ce qu'il comporte : 

- des moyens pour ponderer, par un premier coefficient de 
30 ponderation (agj), au moins un echantillon de rang impair (y2n+m ) voisin 
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d'un echantillon pair en cours de modification, de fagon a obtenir un echantillon 
de rang impair pondere (a 0t j.y2n+mj )• 

- des moyens pour modifier au moins un echantillon de rang pair 
(y 2n ) a partir d'au moins un echantillon de rang impair pondere (<xojy2n+mj ). 

5 - des moyens pour ponderer, par un deuxieme coefficient de 

ponderation (Poj). des echantillons de rang pair (y2n~y2n+2) voisins d'un 
echantillon impair en cours de modification, de fa?on a obtenir des echantillons 
de rang pair ponderes (Poj.(y2n -y2n+2)). et 

- des moyens pour modifier au moins un echantillon de rang impair 
10 (y2n+i) a partir d'au moins un echantillon de rang pair pondere 

(Po,j-(y2n-y2n+2))- 

29. Dispositif de filtrage selon la revendication precedente, 
caracterise en ce que le deuxieme coefficient de ponderation (Poj) est fonction 

du premier coefficient de ponderation (a 0 j). 

15 30. Dispositif de filtrage selon la revendication precedente, 

caracterise en ce que le deuxieme coefficient de ponderation (Po,j) depend du 

premier coefficient de ponderation (aoj) comme suit : 



Poj =nrij/ 



1-2£ao,j 

i=0 j 



ou a 0l i designe le premier coefficient de ponderation, p 0 j designe le 
20 deuxieme coefficient de ponderation, i et j sont des entiers et nij est une valeur 

definie par la recurrence mo =(-1) L ° et rrij = -m^, L 0 etant un entier 
predetermine. 

31. Dispositif de filtrage selon Tune quelconque des revendications 
26 a 30, caracterise en ce que, a chaque iteration, I'echantillon de rang impair 
25 (y2n+mj ) voisin de I'echantillon pair en cours de modification est 

alternativement Techantillon de rang immediatement inferieur (y2n-i) ou 
immediatement superieur (y2n+i). 
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32. Dispositif de filtrage selon Tune quelconque des revendications 
26 a 31, caracterise en ce qu'il comporte des moyens de filtrage 
supplementaires incluant des moyens de ponderation par un troisieme 
coefficient de ponderation (y). 
5 33. Dispositif de filtrage selon la revendication precedente, 

caracterise en ce que le troisieme coefficient de ponderation (y) est fonction 
d'un coefficient de ponderation utilise en amont desdits moyens de filtrage 
supplementaires, comme suit : 

Y = -1/(2P 0i l 0 -i) 

10 ou y designe le troisieme coefficient de ponderation, L 0 est un 

parametre predetermine et Po,i_o-1 d ©signe le coefficient de ponderation utilise 

en amont desdits moyens de filtrage supplementaires. 

34. Dispositif de filtrage selon Tune quelconque des revendications 
26 a 33, caracterise en ce que le signal numerique d'entree (x n ) represente une 

15 image. 

35. Dispositif de filtrage numerique adapte a transformer un ou 
plusieurs signaux numeriques d'entree (y n ) en un signal numerique de sortie 
(x n ), lesdits signaux d'entree comportant des echantillons de rang pair (x 2n ) et 
des echantillons de rang impair (x 2n +i), ledit dispositif de filtrage comportant 

20 - au moins un moyen de ponderation, 

- des moyens pour modifier des echantillons de rang impair (x 2n +i) 
par une fonction d'echantillons de rang pair ponderes (Po,j-( x 2n ~ x 2n+2 )) > 

- des moyens pour modifier des echantillons de rang pair (x 2n ) par 
une fonction (R) d'echantillons de rang impair ponderes (ocoj^n+mj )> 

25 lesdits echantillons ponderes etant fournis par ledit moyen de 

ponderation, lesdits moyens de modification fonctionnant de fagon iterative, de 
fa<pon a modifier au moins une fois des echantillons de rang impair (X2n+i) puis 
au moins une fois des echantillons de rang pair (x 2n ), 

ledit dispositif de filtrage etant caracterise en ce qu'au moins un desdits moyens 
30 de ponderation regoit en entree la difference entre deux echantillons de rang 
pair consecutifs. 
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36. Dispositif de filtrage selon la revendication precedente, 
caracterise en ce que lesdits moyens pour modifier des echantillons de rang 
pair (x 2n ) sont disposes en aval desdits des moyens pour modifier des 
echantillons de rang impair (x 2n +i). 

37. Dispositif de filtrage selon la revendication 35 ou 36, caracterise 
en ce qu'il comporte : 

- des moyens pour ponderer, par un quatrieme coefficient de 
ponderation (Poj), des echantillons de rang pair (x 2n -x 2n+2 ) voisins d'un 

echantillon impair en cours de modification, de fa?on a obtenir des echantillons 
de rang pair ponderes (P 0( j-( x 2n - x 2n+2 )) . 

- des moyens pour modifier au moins un echantillon de rang impair 
(x 2n +i) a partir d'au moins un echantillon de rang pair ponder^ 
(Poj-( x 2n - x 2n+2)). 

- des moyens pour ponderer, par un cinquieme coefficient de 
ponderation (a 0(J ), au moins un echantillon de rang impair (x 2n+m . ) voisin 

d'un echantillon pair en cours de modification, de fagon a obtenir un echantillon 
de rang impair pondere (ocq j.y^+m. ) et 

- des moyens pour modifier au moins un echantillon de rang pair 
(y 2n ) a partir d'au moins un echantillon de rang impair pondere (a 0 j.x 2n +mj )■ 

38. Dispositif de filtrage selon la revendication precedente, 
caracterise en ce que le quatrieme coefficient de ponderation (p 0 j) est fonction 

du cinquieme coefficient de ponderation (a 0 j). 

39. Dispositif de filtrage selon la revendication precedente, 
caracterise en ce que le quatrieme coefficient de ponderation (p 0 j ) depend du 

cinquieme coefficient de ponderation (a 0 j) comme suit : 



Poj = m j 7 



f J ^ 
i=0 
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ou oto.i designe le cinquieme coefficient de ponderation, p 0 j designe 
le quatrieme coefficient de ponderation, i et j sont des entiers et rrij est une 

valeur definie par la recurrence m 0 =(-1) L ° et rrij = -m^, L 0 etant un entier 
predetermine. 

5 40. Dispositif de filtrage selon Tune quelconque des revendications 

35 a 39, caracterise en ce que, a chaque iteration, I'echantiNon de rang impair 
( x 2n+mj ) voisin de I'echantilion pair en cours de modification est 

alternativement Techantillon de rang immediatement inferieur (x 2n -i) ou 
immediatement superieur (x 2 n+i). 
10 41. Dispositif de filtrage selon Tune quelconque des revendications 

35 a 40, caracterise en ce qu'il comporte en outre des moyens de filtrage 
supplementaires incluant des moyens de ponderation par un sixieme coefficient 
de ponderation (y). 

42. Dispositif de filtrage selon la revendication precedente, 
15 caracterise en ce que le sixieme coefficient de ponderation (y) est fonction du 

coefficient de ponderation utilise en aval desdits moyens de filtrage 
supplementaires, comme suit : 

Y = -1/(2p 0i Lo-l) 

ou y designe le sixieme coefficient de ponderation, L 0 est un 
20 parametre predetermine et Po,i_o-i designe le coefficient de ponderation utilise 

en aval desdits moyens de filtrage supplementaires. 

43, Dispositif de filtrage selon Tune quelconque des revendications 
35 a 42, caracterise en ce que le signal numerique de sortie (x n ) represente une 
image. 

25 44. Dispositif de filtrage selon I'une quelconque des revendications 

26 a 43, caracterise en ce que lesdits moyens de modification comportent des 
moyens pour appliquer une fonction d'approximation (R). 

45. Dispositif de filtrage selon la revendication precedente, 
caracterise en ce que la fonction d'approximation (R) est la fonction identite. 



46. Dispositif de filtrage selon la revendication 44, caracterise en ce 
que la fonction d'approximation (R) est une fonction d'une variable reelle qui 
fournit rentier le plus proche de la variable. 

47. Dispositif de filtrage selon la revendication 44, caracterise en ce 
que la fonction d'approximation (R) est une fonction d'une variable reelle qui 
fournit le premier entier inferieur a la variable. 

48. Dispositif de filtrage selon la revendication 44, caracterise en ce 
que la fonction d'approximation (R) est une fonction d'une variable reelle qui 
fournit le premier entier superieur a la variable. 

49. Dispositif de filtrage selon la revendication 44, caracterise en ce 
que la fonction d'approximation (R) est une fonction d'une variable decomposer 
en sous-variables dont la somme est egale a la variable, qui fournit une somme 
de valeurs approchees des sous-variables, chacune des valeurs approchees 
des sous-variables etant, soit une fonction d'une variable reelle qui fournit 
rentier le plus proche de la variable, soit une fonction d'une variable reelle qui 
fournit le premier entier inferieur a la variable, soit une fonction d'une variable 
reelle qui fournit le premier entier superieur a la variable. 

50. Dispositif de traitement de signal (2,5), caracterise en ce qu'il 
comporte un dispositif de filtrage selon Tune quelconque des revendications 26 
a 49. 

51. Dispositif de traitement de signal (2,5) incluant au moins deux 
dispositifs de filtrage selon Tune quelconque des revendications 26 a 49, le 
signal de sortie d'un des dispositifs de filtrage etant le signal d'entree de I'autre 
dispositif de filtrage. 

52. Appareil numerique, caracterise en ce qu'il comporte un dispositif 
de traitement de signal selon I'une quelconque des revendications 25, 50 et 51 . 

53. Appareil photographique numerique, caracterise en ce qu'il 
comporte un dispositif de traitement de signal selon I'une quelconque des 
revendications 25, 50 et 51. 

54. Procede de codage, caracterise en ce qu'il comporte des etapes 
adaptees a mettre en oeuvre un procede de filtrage selon I'une quelconque des 
revendications 1 a 24. 
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55. Dispositif de codage, caracterise en ce qu'il comporte au moins 
un dispositif de filtrage selon Tune quelconque des revendications 26 a 49. 

56. Procede de compression de signaux numeriques, caracterise en 
ce qu'il comporte des etapes adaptees a mettre en oeuvre un procede de 
filtrage selon Tune quelconque des revendications 1 a 24. 

57. Dispositif de compression de signaux numeriques, caracterise en 
ce qu'il comporte au moins un dispositif de filtrage selon Tune quelconque des 
revendications 26 a 49. 
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